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CHAPTER I 


INTRODUCTION 

!.l Combustor Design and Development 

The combustor of the gas turbine engine illustrated in Figure I con- 
tains high intensity combustion and, as far as possible, must burn fuel 
completely. Cause little pressure drop, produce gases of nearly uniform 
temperature, occupy small volume, and maintain stable combustion over a 
wide range of operating conditions (I). Efforts are continually being 
expended to p^’oduce efficient clean combustion. The designer has a for- 
midable problem in aerothermochemi stry , and more accurate analysis pro- 
cedures can provide a route which leads to the accomplishment of design 
objectives more quickly and less expensively than current practice per- 
mits (2, 3). 

In design situations, the engineer has to seek an optimum path be- 
tween alternatives of, for example, efficiency and pollution. The gener- 
al aim of most research investigation is to provide information which is 
useful to designers by "characterizing" or "modeling" certain features 
of the phenomenon in question. Investigations may be theoretical or ex- 
perimental; the two approaches are complementary. Until recently, design 
ers have relied heavily on the experimental approach, but traditional 
design methods (empirical formulas and experimentation) are now being sup 
plemented by analytical methods (numerical solution of the appropriate 
governing partial differential equations). The mathematical modeling 
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approach is now finding favor and is being used to suppiement existirg 
design procedures. Current combustor design and deveiopment problems, 
the needs of the combustion engineer in practice, and proposed research 
tasks which will assist in the attainment of design objectives are becom- 
ing clear. Improvements and new developments (both experimental and the- 
oretical) can and should be made, theoretical modeling oeing aided by 
specific carefully chosen experiments (i») . 

1.2 Practice* Needs and the Present Research 

Many details of combustor phenomena remain poorly understood at the 
macroscopic as well as the microscopic level. Knowledge of these will 
undoubtedly reduce combustor develop>nent time and cost substantial iy, and 
lead to increases in performance. Some of these research needs are: 

1. Physical processes--turbulence, radiation, combustion, and mul- 
t i phase effects. 

2. Computer programs — 0-, I-, 2-, and 3-dimensionai approaches in 
steady state and transient cases. 

3. Unresolved probiems--ef feet of swirl, recirculation, and wail 
proximity on turbulence; turbulence-reaction interaction; and multiphase 
si mu I at ion. 

The present research work is concerned with complementary experimen- 
tal and theoretical studies on a specific subproblem of this general prob- 
lem and is described in a recent paper (S) . The specific problem being 
investigated experimentally is concerned with steady turbulent flow in 
axisymmetric geometries, under low speed and nonreacting conditions — a 
study area highlighted recently as a fundamental research requirement in 
combustion modeling (2, 6, 7). The particular problem is concerned with 


turbulent flow of e given turbulence diitribution in a round pipe enter* 

ing «n expansion into another round pipe, as illustrated in Figure 2. 

The in*coffling flow may possess a swirl component of velocity via passage 

through swirl vanes at angle $ (equal approximately to tan * „) ) • 

and the side-wall may slope at an angle a, to the main flow direction. 

The resulting flowfleld domain may possess a central toroidal recircula* 

tion zone in the middle of the region on the axis, in addition to the 

possibility of a corner recirculation zone near the upper corner provoked 

by the rather sudden enlargement of the cross-sectional area. Of vital 

concern is the characterization of flows of this type in terms of the 

effects of side-wall angle a, swirl vane angle turbulence intensity k, 

in 

of the inlet stream, and expansion ratio D/d on the resulting flowfield 
in terms of its time-mean and turbulence quantities, ^uch problems have 
received little attention, yet there is a definite need for work in this 
area even under nonreacting flow conditions. 

1.3 Theoretical Investigation 

A mathematical solution of the flowfield of interest should provide 
results more economically, quickly, and correctly than possible by other 
means (for example, experiments on real-life systems or nx>dels). In order 
to achieve this, the model should simulate the flow in all Its important 
respects (geometry, boundary condi t ions , phys i cal properties of gases, tur- 
bulence, etc.) and provide a means whereby the governing equations may be 
solved. Mathematical models of steadily increasing realism and refine- 
ment are now being developed, both in the dimensionality of the model 
(together with the computational procedures) and in problems associated 
with the simulation of the physical processes occurring. Clearly there 


3 




•r« two areas of difficulty: lha simulation and the solution. Recent 

publications discuss these difficulties in the context of combustor flow- 
field predictions (k, Si 8, 9i 10). 

An advanced prediction computer program, which solves the governing 
partial differential equations in finite difference form, has been devel- 
oped. A detailed discussion of this program, including a user-oriented 
section, has been prepared recently (II). Therefore, extensive discus- 
sion of this program in Chapter 1 1 1 of this docu«nent Is not necessary and 
only highlights of the technique are included. 

I.A Experimental Investigation 

Prediction work in the form of computing the fiowfield is being com- 
plemented by an associated experimental study which Is used to evaluate 
the final predictive capability. The experiments, under way at Oklahoma 
State University, are concerned with measuring the effects of swirl and 
side-wall angle on streamlines as well as mean flow and turbulence para- 
meters in nonreacting flow. Use is made of the confined Jet facility 
shown schematically in Figure 3« 

The facility has an axial flow fan whose speed can be changed by 
means of a varidrive mechanism. Numerous fine screens and straws produce 
flew in the settling chamber of relatively low turbulence intensity. The 
contraction section leading to the test section has been designed to pro- 
duce a minimum adverse pressure gradient on the boundary layer and thus 
avoid unsteady problems associated with local separation regions. The 
sudden expansion consists of a 15 cm diameter circular Jet nozzle, exit- 
ing abruptly into a 30 cm diameter test section, as shown in Figure 3. 
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Th« substantial size of this test model provides excellent probe resoiu* 
tion for five-hole pitot measurements. The test section is constructeo 
of plexiglass to facilitate flow visual izat Ion. The side-wall angle 
and swirl vane angle i are variable. The side-wall angle is set by in- 
serting one of three blocks with a side-wall angle a of 90, 70, or hS de- 
grees. The swl<*l vane angle Is continuously variable with discrete val- 
ues of ^ ■ 0, h 9, and 70 degrees being chosen for the present study. 

The experinients are concerned with nonreacting flow, and the following 
two aspects of the study are addressed in the present thesis: 

1. Photography of neutral I y-buoyant helium-filled soap bubbles, 
tufts, and smoke, so as to characterize the time-mean streamlines, recir- 
culation zones, and regions of highly turbulent flow. 

2. Five-hole pitot probe pressure measurements to determine time- 
mean velocities u, v, and w. 

Complete details of these experimental studies appear in Chapter IV. 

I. 5 Swirl Flow Concepts 

Host practical combustion equipment includes the phenomenon of swirl 
and It is %'nvenient to recall here some of the basic concepts of its gen- 
eration and characterization (12). Swirling flows result from the appli- 
cation of a spiraling (notion, with a swirl velocity component (also known 
as a tangential or azimuthal velocity component) being imparted to the 
flow via the use of swirl vanes, in an axi al -pi us-tangent iai entry swirl 
generator or by direct tangential entry Into the cofnbustion chamber. Ex- 
perirnental studies show that swirl has large scale effects on flowfields; 
jet growth, entr-*nment, and decay (for Inert jets) and flame size, shape, 
stability, and combustion intensity (for reacting fUiws) are affected by 
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th« degree of swirl imparted to the flow. This degree of swirl usually 
Is character i zed by the swirl number S, which is a nondimensional number 
representing axial flux of swirl momentum divided by axial flux of axial 
momentum times equivalent nozzle radius. That is 


S 


“e 

ir~m 


(i.i) 


where 


Gg - (puw ♦ pu'w' ) r dr 

is the axial flux of swirl momentum, including the x6-direction 
turbulent shear stress term; 

2 2 

(pu pu' ♦ (p - p^)) r dr 


is the axial flux of axial momentum, including tne x-direction 
turbulent normal stress term and a pressure term; d/2 is the nozzle 
radius; and u,v,w are velocity components in (x,r,8) cylindrical 
polar coordinate directions. 

In a free Jet in stagnant surroundings, G and G. are constants, that is, 

X o 

invariants of the jet, as may be deduced from radial integration of the 
appropriate governing Reynolds equation (13)* 

It Is especially convenient to be able to relate the angle of a 
swirl vane pack to the swirl number produced. In this context, for com- 
parison purposes it is worth noting that swirl vane angle ^ and swirl 
number S are related approximately by 


S 


2 

^ \i - (d^/d)^y 


tan^) 


( 1 . 2 ) 


« 
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wh«r« d and d. are noazia and vana pack hub diimatara, ratpactivcly. 1nl« 
n 

ralatlonthip folloMS from aaaumptlona of plug flow axial valocity u^ in 

the annular region, and very thin vanes at constant angle # to the main 

direction to imparting a constant swirl velocity " u tanf to the flow. 

o o 

From the and definitions, Integration 0 '.*er the range (■ to 

R (■ d/2) gives 



G 


X 




(1.3a) 

(1.3b) 


from which the quoted relationship follows, in the case of a hubless 
swirler (or one for which d./d is very small) the expression simplifies 


S ■ j tan^ ( 1 .4) 

so that vane .ingles of 15, 30, 45, bO, 70, and 80 degrees, for examole, 
correspond to S values of approximately 0.2, 0.4, 0.7, 1.2, 2.0, and 4.0. 
Hera 100 percent efficiency also is assumed for the swirl vanes, which 
deteriorates as the vane angle increases, in fact, flat blades are not 
perfect, and indeed the assumption of flat axial and swirl velocity radial 
profiles invTNidiately after the swirler is not correct, as may be seen frCNr< 
experimental measurements in Chapter VI of this thesis and elsewhere (12, 
13. 14). 


1.6 The Present Contribution 

The basic objectives of the present research include the development 
of a computer code for the computation of swirling, ax I symmetric, confined 
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J«t floMft ufting computational procadurat and tha k>c 

turbulanca modal. Furthar objactivat antall mean floMflald maaturamantft 
and floM vUual Izat Ion to charactarlia tbit particular flow and to avalu* 
ata tha davciopad pradictlon procadura. Tha computations and experiments 
hava bean performed for six basic confined Jat flowflalds, all of which 
conform to tha schematic diagram of Figure 2 and hava an o«>^;rsnslon ratio 
0/d ■ 2. Tha parameter combinations which define these fiowflaids are 
side-wall angle o ■ 90 and kS degrees and swirl vane angle g ■ 0, kS, and 
70 degrees. The computational techniques are discussed in Chapter III. 

The experimental program consisted of performing mean velocity mea- 
surements wi th a f I ve-hole pi tot probe and flow visualization ofneutrally- 
buoyant soap bubbles, tufts, and smoke. Details of the experimental 
methods ere presented in Chapter IV. The pitot probe measurements re- 
sulted in maps of axial, radial, and swirl velocity componerits while the 
flow visual Izat ion provided streamline patterns and recirculation z^.e 
geometries. Results of the overall flowfielc* characteristics highlighted 
by the recirculation region geometries, both predicted and visualized, 
are presented in Chapter V. Chapter VI contains results of predictions 
and measurements of tlme-meen velocity coi'>i>onents . Finally, Chapter VII 
summarizes the conclusions of the present investigation and gives recom- 
mendations for further work. 



CHAPTER I I 


REVIEW OF PREVIOUS AXISYHMETRIC 
RECIRCULATING FLOW STUDIES 

Previous experimental and theoretical work In axisymmetric recircu- 
lating flows provides an Important background for this Investigation. 
Information regarding both swirling and nonswirling, nonreacting, turbu- 
lent flows has been extensively surveyed. Noteworthy results are Includ- 
ed In the following summary where the emphasis is on nonreacting turbu- 
lent swirling flows. 


2.1 Experimental Work In Nonswirling Flows 

The effect of expansion angle a was investigated by Chaturvcdi (15) 
who measured mean and turbulent flow quantities behind expansion half- 
angles o' 15. 30, 45, and 90 degrees for air flow in a pipe. The expan- 
sion diameter ratio was 2.0. Measurements of velocity head in regions of 
high turbulence Intensity and where the direction of the velocity vector 
was unknown were n.ade with a 2.5 mm diameter pitot tuba. Mean velocity 
was also measured with a constant temperature hot-wire anemometer using a 
single wire. A cross-wire was used to measure all the Reynolds stresses. 

Krall and Sparrow (16) utilized an orifice in an electrically heated 
tube to create flow separation. It was presumed that the peak Nusselt 
number occurred at the reattachment point. Their results agreed with 
those of Phaneuf and Netzer (17) In that reattachment length was unaffected 
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by Reynolds number in the turbulent range. Also, they found that peak 
Nusselt number became spread out Into a broad zone for strong separations 
caussd by large step heights. This transition from a point to a broad 
zone occurred between h/D ■ 0.16 and h/D ■ 0.25, where D is the pipe dia* 
meter and h is the step height. Also, It was found that Increasing the 
step height moved the peak point slightly downstream. 

Back and Roschke (l8, 19) investigated the reattachment iength be* 
hind a pipe step by visual observation In water flows. The diameter ex* 
pension ratio was 2.6 and a conical contraction section Just upstream 
reduced the boundary layer thickness. Small 0.76 mm diameter holes were 
located along the larger diameter tube one step height apart. Reattach* 
ment locations were determined by slowly metering dye through these holes 
and observing whether it moved upstream or downstream. 

Wall static pressure was measured by Phaneuf and Netzer (17) as well 
as mean axial velocity behind abrupt pipe steps wherein the step heights 
were 0.303 and 0.199 outlet diameters. They found that the reattachri-ient 
zone spreads out with both increasing Reynolds number (for a given step 
height) and increasing step height provided the Reynolds number is in the 
turbulence transition range. It was further concluded that the point of 
maximum heat transfer does not coincide with the zone of reattachment, 
and furthermore, it is dependent upon the inlet mass flux. 

Ha Mlnh and Chassaing (20) studied the restructuration of a turbu- 
lent pipe flow past three types of geometric perturbations, one of which 
was the abrupt pipe expansion. Agreement between hot-wire and pitot 
probe measurements of time-mean axial velocity was exhibited. Owing to 
probe resolution considerations, only normal or inclined single wire 
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probe» w«re used. All Reynolds stresses were measured, and the rotated 
inclined hot'wire technique was utilized. 

Both a one*color and two-color I aser-Ooppler system were utilized by 
Hoon and Rudinger (21) to measure mean axial velocities behind an abrupt 
pipe step. The diameter expansion ratio 0/d was l.i«3* Velocity reverses 
in the recirculation zone were easily detectable through use of a Bragg 
cell to shift the frequency of one probe beam relative to the other. 

An experimental technique was developed by Johnson and Bennett (22) 
for measuring mean and fluctuating velocity and concentration distribu- 
tions along with their cross correlations. Measu>'ements were obtained 
using a laser-Dcppler velocimeter to measure velocities, and a laser- 
induced fluorescence technique to measure the concentration of fluorescent 
tracer particles. Such measurements were displayed for coaxial Jets dis- 
charging into a larger axisymmetric duct with D/d ■ 2.07. 

2.2 Experimental Work in Swirling FlowV' 

Both an aircraft and an industrial type combustor were investigated 
by Hiett and Powell (23) under reacting and nonreacting conditions. A 
water-cooled three-dimensional pitot probe was developed. Mean velocity 
profiles, recirculation boundaries, and temperature contours were present- 
ed. Alio, cotton tufts were employed to visualize the nonreacting flow. 

Extensive studies of free and enclosed swirling jets from vane swirl- 
ers were undertaken by Mathur and MacCallum (24). In the enclosed case, 
the expansion ratio was approximately D/d ■ 2.5 into a square-sectioned 
test section. An evaluation of swirler efficiency was conducted for both 
annular (with hub) and hubless swirlers. Extensive presentation is given 
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of velocity profiles and recirculation zone sizes as functions of the de* 
gree of inlet swi rl . 

Laser velocimeter msasurements were made by Owen (25) In the Initial 
mixing region of a confined turbulent diffusion flame burner. The experi' 
mental facility consisted of an axisymmetric combustor in which a central 
gaseous fuel stream mixed with a high speed coaxial annular air stream. 

The expansion diameter ratio was approximately 1.3 end swirl numbers of 
0, 0.3i end 0.6 were chosen to contrast swirling with nonswirling combust* 
ing flowfields. Measurements were presented for axial and tangential mean 
velocity profiles along with the rms and probability density distributions 
of velocity fluctuations. 

Three mean velocity components and normal stresses weie measured with 
a laser anemometer by Baker et al. (26) in an experimenia! facility of a 
small scale axisymmetric furnace. The test chamber has an expansion dia- 
meter ratio of 3.33 and contains a constricted outlet. Measurements were 
obtained in both an isothermal air flow and a combusting mixture of natu- 
ral gas and air, with exit swirl numbers of 0 and 0.52. The boundary and 
field measurements of three components of time-mean velocity and corre- 
sponding normal stresses were designed to be of good use In the evaluation 
of turbulent flow prediction procedures. It was clear from the results 
that regions of recirculation under reacting flow conditions differed sub- 
stantially from those under nonreacting conditions. Surprisingly, they 
were found to be in general larger under reacting flow conditions. The 
turbulence was found to be far from isotropic over most of the flow field. 

Beltaqui and MacCailum (27> 28) studied premixed vane swirled flames 
in furnaces with expansion ratios D/d ■ 2.5 and 5.0, where U and d are 
chamber and nozzle diameters, respectively. They found that a given vane 
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iwirler can produce diff*;rent velocity profiles as D is varied, and that 
Che size and shape of the central recirculation zone are primarily func- 
tions of D rather than d. A modified swirl number using D rather than d 
was suggested as being more characteristic. 

Mean velocity measurements were obtained by Vu and Gouldin (29) for 
coaxial jet mixing in a pipe under coswirl and counterswirl conditions. 
Both the central and annular flows had a swirl number of approximately 
O.S. In the counterswirl case, a recirculation zone formed on the axis 
with very steep velocity gradients. In the coswirl case, no recircula- 
tion was observed. 

Syred et al . (30) concentrated on turbulence measurements in strong- 
ly swirling flows via a single-wire six-orientation hot-wire technique. 
Local turbulence intensities were found to be extremeiy high in and near 
the central recirculation zone. Measurements of all six turbulent stress 
components showed strong variations of absolute turbulent kinetic energy 
ievels and strong nonisotropy of the stresses and associated turbulent 
viscos i ty. 

Further work of less importance to the present study is readily 
available (31. 32, 33. 3^*). 

2.3 Theoretical Work in Nonswirling Flows 

Laminar flow over a pipe step has been computed for the diameter ex- 
pansion ratio 2.0 by Macagno and Hung (35). Computations for Reynolds 
numbers up to 200 were undertaken using stream function and vorticity 
variables. Flow visualization results were in good agreement with pre- 
dicted values. 
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Runcha) at al. (36) also predicted the flow downstream of a pipe 
step, solving transport equations for stream function, vorticity, enthal- 
py, and turbulence energy. The turbulence length scale was empirically 
determined to ensure that the reattachment point was accurately predicted. 
The effects of the outlet boundary condition are significant for only a 
short distance upstream, for example, one pipe diameter. The Stanton num- 
ber was correctly predicted to vary as the -O .38 power of the Reynolds 
number. 

As previously mentioned, Ha Minh ard Chassaing (20) investigated the 
flow behind a pipe step. They obtained good agreement between measured 
and predicted values considering the extent of the flow perturbation. 

The k-c turbulence model was employed for most of the computations, al- 
though the five-equation Reynolds stress transport model predicted the 
centerline kinetic energy more accurately. 

Other predictive studies on nonswirling flow are discussed in sec- 
tion 2. A, where the emphasis Is on swirling flows and in which the non- 
swirling case appears as a degenerate case. Selected theoretical work on 
turbulence modeling is also useful (37, 38, 39, ^0). 

2.k Theoretical Work in Swirling Flows 

Roberts (41) investigated similar axisymmetric, swirling flowflelds 
by developing a finite difference code incorporating the k-W turbulence 
model while using stream function and vorticity variables. Poor perform- 
ance of the k-W model led to the crude approach of specifying constant 
effective viscosities where ■ “lO. Also, static wall 

pressures were measured and visualization information recorded for a 
water model experiment which gave general agreement with predictions. 
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Khalil et al. (42) predicted the flowfleld of several two*dlmenslonal 
model furnaces for which various measured quantities are available using a 
similar finite difference algorithm. As detailed knowledge of inlet boun* 
dary conditions are rarely known, a parametric study was undertaken, and 
results indicated that the inlet velocity profile has negligible effect. 

Also, provided that total enthalpy of the inlet flow is held constant, 
changes in the radial temperature profile exhibited a very small effect. 

However, it was found that inlet distributions of k and c exert an appre- 
ciable Influence, and of the three inlet k profiles tested for the flow- 

2 

field measured by Baker et al. (26), values for k/u of 0.009 for near 

o 

wall nodes and 0.003 for central nodes gave closest agreement for axial 
velocity. Nonswirling flows are more sensitive to inlet profiles and 
thus the above numerical tests were conducted without swirl. Also, no 
significant differences could be observed for various outlet boundary con- 
ditions tested for axial velocity, enthalpy, mixture fraction, and concen- 
tration fluctuation. 

An axisymmet r i c , combusting, swirling flow computer code was used by 
Hutchinson et al. (43) to compare predicted and measured quantities in 
the furnace test facility studied by Baker et al. (26) except at higher 
flow rates. The furnace consisted of a coaxial burner with swirling annu- 
lar air flow and a constricted exit. Measurements of mean velocity and 
normal stress in both axial and tangential directions were generally pre- 
dicted along with mean temperature and wall heat flux. Similar calcula- 
tions using the same physical information but with a different numerical 
scheme showed that differences of up to 15 percent can be introduced by 
the numerical arrangement, whereas the rms differences were less than 5 ( 

percent. Discrepancies in axial velocity ranged up to 40 percent of the 

15 
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maxitnum In the appropriate plane. The largest discrepancies were in the 
burner wake and It was suggested that the two*equatlon turbulence model 
is inadequate there. 

A finite difference code was developed by Li I ley (44) for inert and 
reacting recirculating flows with strong swirl. The effect of swirl num- 
ber on mean axial and swirl velocities was presented along with predicted 
streamlines in an axisymmetrlc combustor. Also, the effect of swirl vane 
angle on central recirculation zone length and flame length was exhibited. 

Serag-Eldin and Spalding (45) predicted the nonaxi symmetr i c, swirl- 
ing, combusting flow in a can combustor with air entering through a 45*’ 
degree vane angle swirler fn an annular tube enclosing a central jut of 
fuel. Temperature contours were measured with a bare 13 percent Rh-Pt 
versus Pt thermocouple wire of 3*8 x 10 ^ mm diameter. Predicted results 
gave the same trend as measured ones, but quantitative agreement was not 
always acceptable. 

Swirling, combusting flow in axisymmetrlc prechamber dump combustors 
was predicted by Novick et al. (46) using a finite difference procedure. 
Application of the computer code *o specific combustor designs is stress- 
ed, and consideration is given to the relative effect of certain impor- 
tant design parameters: main chamber to prechamber expansion ratio, 

swirler dimensions, double concentric reverse swirler, constricted out- 
let, and fuel concentration profile. The standard base configuration con- 
sists of a 70-degree swirl vane angle inlet flow (stoichiometric mixture) 
which proceeds over a central hub, enters a prechamber, followed by a sud- 
den expansion. Computed streamlines and temperature contours for this 
case are shown. 



16 


OetaiU of the above model Including physical assumptions, governing 
equations, finite difference equations, and solution procedure are found 
In a recent paper by Novick et al. (k7)« The predictive capability Is 
demonstrated for a flowfield with coaxial Jets entering a sudden expan- 
sion. Computed streamline patterns are presented for various *.wirl 
strengths with and without combustion in a prechamber dump combustor. 

Other advancas to prediction procedures have been made elsewhere and 
applied to other two-dimensional axisymmetric problems involving swirl 
flows in expansions of interest to the combustor designer. Habib and 
Whitelaw (k8) compared hot-wire and laser measurements with predictions 
in confined coaxial Jets with and without swirl. 

Predictions of the swirl mixing pipe flow data of Vu and Gouldin 
( 29 ) were made by Srlnivasan and Mongia (k9) for both coswirl and counter- 
swirl cases. The measured inlet profiles were utilized as inlet boundary 
conditions with both the original k-e turbulence model and the version ex- 
tended to include Richardson number effects. The extended model gave bet- 
ter agreement with measurements, but the recirculation zone was somewhat 
elongated. 

Sturgess and colleagues have predicted a variety of flow types with 
and without combustion as part of a continuing validation study of turbu- 
lence and reaction models for aircraft gas turbine combustor application 
(SO, SI, S2). Associated sudden expansion experiments are yielding much 
needed data for this evaluation (22). 

Li I ley et al. (9) discussed the predicted effects of swirl vane angle 
^ and inlet side-wall expansion angle a on axisymmetric idealized gas tur- 
bine combustor flowfields. Also, the predicted effects of a and ^ on both 
the corner and central recirculation zone lengths were presented, under 
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loth reacting and nonraacting flow conditions, so extending earlier iso’ 
thermal studies. 

Other papers presented at a recent conference address turbulent re- 
acting fiows in practical combust io" systems (S3). The computation of 
flows with complex cnemical kinetics problems, complex fluid mechanic 
problems, and complex simulation problems are all considered. 

Fully three-dimensional flowfields are not specifically addressed in 
the present thesis. However, by way of introduction, interested readers 
are referred to the work of Mongla and Reynolds (Sk) , Serag-Eldin and 
Spalding (kS) , and Swithenbank et al. (SS). Mongia and Reynolds (Sk) 
developed and applied an advanced fully three-dimensional prediction tech- 
nique to compute ihe reacting flowfield as part of an analytical desi^r. 
procedure for gas turbine combustion systems. Complexities such as the 
eddy-breakup chemical reaction model, and fuel spray droplet phenomena 
were included. Good qualitative agreement was found to be quite reward- 
ing in developing a number of relatively complex combustion systems. 
Gaseous phase diffusion flames In can-type combustors have been predicted 
(kS). Patterns of velocity and temperature show that qualitative trends 
may be simulated. A three-dimensional two-phase mathematical model of a 
Lycoming combustor has been amalgamated by Swithenbank et al. (55). Some 
30 partial differential equations were solved and encouraging results 
were obtained. 
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CHAPTER I I I 


COMPUTATIONAL TECHNIQUES 

The economical design and operation of practicol fluid firr./ equip- 
ment can be greatly facilitated by the availability of prior predictions 
of the flowfield. These may be obtained by use of a mathematical model 
incorporating a numerical finite difference solution procedure, and sev- 
eral textbooks discuss the approach at length (56, 57i 58). Several re- 
cant review papers also provide useful details (Ik, 37, k2 , kk, kS, k7, 52, 
59). These theoretical studies complement associated experimental stud- 
ies, and their improvement and use can significantly increase understand- 
ing and reduce the time and cost of development work. In the context of 
the present study, the main goals Include improved flowfield simulation 
and turbulence model development, and computer code development plays a 
significant role in this endeavor. 

The computer program STARPIC (acronym for girling t^urbulent axlsym- 
metric recirculating flow in practical jjsothermal £ombustor geometries) 
has been developed and employed in the theoretical aspects of this inves- 
tigation. It is used in a parametric study to predict mean velocity pro- 
files and streamline patterns under various flow conditions. The develop- 
ment of STARPIC began with the nonswirling TEACH (peaching e^lliptlc axi- 
symmetric £haracter istics heurist ical ly) computer program (60) as opposed 
to that of LI I ley (kk) because the latter contains extensive combustion 
calculations which are inappropriate under the present objectives. 
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B«sld«» th« Incorporation of all awirl momantum af facts , numarous special 
faatiiras wara Includad, such as: (I) a panarallzad "stairstap" simulation 
of tha sloping combustor wall, (2) an advancad nonuniform grid systam, 

(3) momentum wall functions derived from a recant axparlmant, (k) use of 
advanced estimates of Initial field values, and (5) streamline calcula* 
tion and plotting. 

Discussion of this finite difference computer program is minimized 
hare because it Is extensively described in a recent NASA report (II). 

The report deals with the complete computational problem, showing how the 
mathematical basis and computat lonal scheme may be translated into an 
effective computer program. A flow chart, FORTRAN IV 1 1st Ing, notes about 
various subroutines, and a user's guide are also supplied. 

3.1 The Governing Equations 

The turbulent Reynolds equations for conservation of mass, momentum 
(in X, r, and 6 directions), turbulence energy k, and turbulence dissipa* 
tion rate r. which govern the two-dimensional, axisymnetric, swirling, 
steady flow may be taken as previously (^40, k2, kk, k7, ^8, 49, 61). The 
transport equations are all similar and contain terms for convection and 
diffusion (via turbulent flux terms) of a general dependent variable d* 

The source term contains terms describing the generation (creation) 
and consumption (dissipation) of d. Introducing turbulent exchange co- 
efficients and the usual turbulent diffusion-flux (stress - rate of strain • 
type) laws, it can be shown that the similarity between the differential 
equations and their diffusion relations aliows them ali to be put in the 
common form: 
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i <pur») * i (pvrp) - i (rr^ |i) 




(3.U 


Th« forms of the source term are given In Table I and other quant I* 
ties are defined as follows: 


? / ^ ^ / _ Sv \ 
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(3.2a) 

(3.2b) 

(3.2c) 
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(3. 2(1) 


implicit here is the use of the two-equation k-c turbulence model 
(62, 63) and constant Prandtl-Schmidt numbers from which exchange coeffi' 
cients are calculated: 


U ■ C pk /c ♦ w* 

u * 




(3.3) 


(3.M 


The total equation set must be solved for the time-niean pressure p, 
velocity components u, v, and w, etc. Then other useful designer inform- 
ation, like streamline plots, breakaway and reattachment points, recircu- 
lation zones and stagnation points, for example, may be readily deduced. 

3.2 The Finite Difference Formulation 


Solution may be via the stream funct ion-vort icity, or primitive 
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pr«s$ur«-v«loci ty Approach. Whereat the former approach, utao in tht; 

1968 computer program from Imperial Col 'age for example (66), reduces by 
one the number of equations to he solved /..'d eliminates the troublesome 
pressure (at the expense of trouble with the vo'tlcUy equation), the 
preferred approach now Is SIHPLL (mnemonic for Mml •j_mpl Id t method for 
pressure J[l iked e<)uatlons) which focuses attention directly on the let* 
ter variables. 

Because It possesses many advantages, the present work has been 
developed using this new technique, which was Incorporated In the TEACH 
program (60). 

The finite difference equations are solved on a complex mesh system, 
shown in Figure 6. The intersections, the point P for example, of the 
solid lines mark the grid nodes where all variables except the u and v 
velocity components are stored. The latter are stored at points which 
are denoted by arrows (and labeled w and s, respectively) located midway 
between the grid intersections. Details of the special merits of this 
staggered grid system have been reported previously (60, 6S). The dif- 
ferent control volumes C, U, and V, which are appropriate for the P, w, 
and s locations, respectively, are also given in the figure. Thus the 
u- and v-cell bou.idarles lie either exactly on or exactly halfway be- 
tween the grid intersection points P. This improves the accuracy of cal- 
culating pressure gradient effects on the u and v velocities (as pressure 
values now lie directly on appropriate cell faces) and of calculating 
certain ceil boundary fluxes. The flowfleld domain boundaries are posi- 
tioned midway between grid points so that tliey coincide with normal velo- 
cities as in Figure 5. which shows an example of grIJ specification for 
the geometry under consideration. 
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In obtaining the finite difference equations using the primitive 
approach, it is common to use a hybrid formulation, which is a combina- 
tion of the central and upwind differencing scheme, it is applied to 
the rnpresentat ion of the convection and diffusion tenns and provides 
(when needed) the required degree of upstream differencing. In this man- 
ner the finite difference equation for a typical variable becomes 

•* 

where 


a 

P 


I 


j 




sum over 


P 

N, 


S, E, and W neighbors 


for each variable at each point P of a grid system covering the flow 
doma i n . 


3.3 The Solution Procedure 

The flowfield is covered with a nonuniform rectangular grid system; 
typically the boundary of the solution domain falls halfway between its 
Immediate nearby parallel gridlines; and clearly specification of the x 
and r coordinates of the gridlines, together with information concerned 
with the position of the boundary, is sufficient to determine the flow- 
field of interest. The curved or sloping boundaries are simulated by 
means of a stairstep approach with sloping boundary segments where appro 
priate. 

Finite difference equations like Equation (3.5) hold for t 11 vari- 
ables at all internal points of the mesh system, except those adjacent 
to wall boundaries. At these points, the correct boundary influence is 
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Inserted by a linearized source technique, using wall functions so as to | 

avoid detailed calculations in near-wall regions. The effect of swirl on 
wall function specification is handled as follows. The previous ideas 
are extended to find total tangential wall shear stress near boundaries 
(involving x and 6 directions for *n r ■ constant wall, and r and 6 direc- 
tions for an X ■ constant wall). Then appropriate components are deduced 
directly (for the u and w velocities which are tangential to an r ■ con- 
stant wall, and v and w velocities which are tangential to an x* constant 
wail). The effects on u, v, and w momentum equations are incorporated 
via the usual linearized source technique (II). 

The finite difference equations and boundary conditions constitute a 
system of strongly-coupled simultaneous algebraic equations. Though they 
appear linear they are not, since the coefficients and source terms are 
themselves functions of some of the variables; and the velocity equations 
are strongly linked through the pressure. In the solution procedure alge- 
braic equations like Equation (3>5) are solved many times, coefficient 
and source updating being carried out prior to each occasion. The prac- 
tice used here »s to make use of the well-known i^ri -^i agonal matrix algo- 
rithm (ThMA) , whereby a set of equations, each with exactly three unknowns 
in a particular order except the first and last which have exactly two un- 
knowns, may be solved sequentially. In the two-dimensional problem one 
considers the values at gridpoints along a vertical gridline to be unknown j 

i 

(values at P, N, and S for each point P) , and the most recent values at . ' 

each E and W neighbor are considered known. The TuMA Is then applied to 
this vertical gridline. In this manner one can traverse along all lines 
in the vertical direction sequential ly from left to right of the integration 
domain, using a certain degree of underre 1 axat ion to enhance convergence. 
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CHAPTER IV 


EXPERIMENTAL FACILITY AND TECHNIQUES 
4.1 Wind Tunnel 

The Oklahoma State University swirling confined Jet test facility is 
used in the experimental study. A schematic of the overall facility is 
shown in Figure 3> Ambient air enters the low'speed wind tunnel through 
an air filter consisting of foam rubber. Then it flows through an axial 
flow fan whose speed is continuously variable over a wide range in that 
it is driven by a S H.p. varidrive motor manufactured by U.S. Electrical 
Motors, Inc. Next, the flow gradually expands with the tunnel cross sec- 
tion without separation because numerous fine mesh screens are encounter- 
ed along the way. 

The flow conditioning section is next. It consists of: a perforat- 

ed aluminum plate (2 mm diameter holes), followed by a fine mesh screen, 
a 12.7 cm length of packed straws, and five more fine mesh screens. Next 
is the contoured nozzle leading to the test section. This axisymmetric 
nozzle was designed by the method of Morel (66) to produce a minimum ad- 
verse pressure gradient on the boundary layer to av id the flow unsteadi- 
ness phenomenon associated with local separation ms. The area ratio 

of the flow conditioner cross section to that of the nozzle throat is 
approximately 22. 3> Further details are available elsewhere regarding 
the design and construction of the nozzle, flow conditioner, and ideal- 
ized test model which is described in the next section (67, 68). 
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k.2 Test Section 


The test section consists of a swirl vane assembly and an idealized 
combustion chamber model which is mounted onto the throat of the wind 
tunnel nozzle. A photograph of the general arrangement in the laboratory 
is given in Figure 6. The swirl vjne assembly contains ten brass vanes 
which are individuaiiy adjus'.aole for any vane angle The diameter of 
the flow passage within the swirler matches that of the contour nozzle 
throat which is approximately 15 cm. A solid hub of ^ cm diameter is 
located at the center of the swirler with a streamlined nose facing up- 
stream, as seen in Figure 7. The downstream end is simply a flat face, 
simulating the geometric shape of a typical fuel spray nozzle. Figure 8 
shows the swirler assembly from the downstream side. A special Jig, also 
seen In the figure, is employed to set each vane at the desired angle. 

The idealized combustion chamber model is composed of an expansion 
block and a long plexiglass tube. The expansion block is a 30 cm diameter 
disk of oak wood which is mounted to the downstream section of the swirl 
generator. The block has a 15 cm diameter hole centered on its axis 
through which the air entering the model flows smoothly. The downstream 
face of the expansion block has been shaped to provide the desired flow 
expansion angle a which is shown in Figure 9. There are currently three 
interchangeable expansion blocks and the appropriate choice gives a ■ 90, 
70, or 45 degrees. Figure 9 shows a view of the swirler with two of the 
expansion blocks. 

The actual combustor of a gas turbine engine, which is illustrated 
in Figure 1, is idealized in the present study as there are no film cool- 
ing holes or dilution air holes, and the simulated chamber wall consists 
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of a straight pip«. The inside diameter of the pipe is 30 cm and the 
iength is approximateiy I2S cm. The substantiai size of the test model 
provides excellent probe resolution for five-hole pitot probe measure- 
ments, and it also enhances flow visualization. Plexiglass was the 
chosen material for the test section pipe In order to facilitate flow 
visualization. It is supported by a wooden cradle which is mounted on a 
table. 

The test section is carefully aligned with the wind tunnel so that 
the test section and wind tunnel centerlines are col inear. The swirl vane 
assembly and the expansion block are always securely centered on the wind 
tunnel centerline, as they are both positioned on the discharge end of 
the contoured contraction with dowel pins. The test section inlet is 
centered since it slips snugly over the expansion block which was care- 
fully machined with a lathe. The test section discharge, however, is 
aligned with a low-power Spectra-Phys i cs StabiUte, model 120, laser beam 
which passes through its center. A plexiglass disk which fits tightly in 
the model discharge identifies this center. 

A removable length scale exhibiting axial distances downstream is 
mounted above the test section for flow visualization photography. 

4.3 Flow Visualization 

Flow visualization experiments using three different techniques were 
conducted during this investigation and each is described in the follow- 
ing subsections. Such experiments serve as an important investigative 
tool by immediately providing a clear view of the flowfield. During the 
initial phase of this study, it was employed: (1) to examine the overall 

flowfield; (2) to detect unforeseen difficulties with measurements; (3) 
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to check the gross character istics of the coniputer predictions; (4) to 
locate possible regions of interest in which additional probe measure* 
inents might be desired; and (5) to help interpret measured data. The un- 
foreseen difficulties with measurements mentioned above include time- 
dependent instability effects, regions of severe turbulence, or probe 
i nterference. 

Probably the most widely used flow tracer to date is smoke. Smoke 
particles are very small, typically less than 1.0 x 10 ^ mm in diameter. 
Hence they follow air motion quite well; however, they cannot be viewed 
or photographed individually. Thus one can only create a relatively 
dense cloud of particles and view the motion of its edges. The cloud 
will diffuse if the air flow is turbulent, the edges becoming indistinct 
and then disappearing as the particle concentration drops. An attractive 
flow visualization technique for both detailed as well as overall charac- 
terization in strongly turbulent flows employs neutral ly-buoyant helium- 
filled soap bubbles. Although much larger than particles of smoke, soap 
bubbles are almost ideal because they trace particle paths and have the 
average density of air (5). 

4 . 3.1 Neutral ly-Buoyant Soap Bubbles 

Apparently soap bubbles were first used as flow tracers by H. C. H. 
Townend, whose work was reported by Lock (69) in 1928. Coal gas was used 
to make the bubbles neut ra 1 ly-buoyant for use in a propeller flowfield 
study. Redon and Vinsonneau (70) injected air and also hydrogen into 
soap bubbles in a wind tunnel and concluded that air-filled bubbles gave 
essentially correct trajectories. Some statistical measurements of tur- 
bulence using bubbles in a wind tunnel were reported by Kampe de Feriet 
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(70. Of all the previous work with soap bubbles, that nost related to 
the present study was conducted by Owen {J2) , who employed them to study 
the flowfield in a cylindrical vortex tube. 

The bubble flow visualization system used in this study is composed 
of (I) the bubble generator, (2) the light source, and (3) two cameras 
and tripods. The arrangement of some of these components is schematical- 
ly illustrated in Figure 10. The bubble generator is model 3, manufactur- 
ed by Sage Action, Inc. It consists of the Injector, shown in Figure II, 
in which the bubbles are actually formed, and a console that supplies the 
constituents to the injector. Only the injector is situated in the tun- 
nel. It is located far enough upstream of the test section to avoid any 
disturbance to the flow pattern under study. The injector currently em- 
ployed is 1.6 cm in diameter and approximately 11.5 cm long with a I cm 
diameter stem perpendicular to the main body. Compressed air at a gauge 
pressure of 4.1 bar is supplied to the injector through this stem, while ' 

helium at a gauge pressure of 1.4 bar and bubble film solution are intro- 
duced through smaller tubes directly upstream of the stem. The bubble 
film solution is a special soap mixture, which is commercially available 
from Sage Action, Inc. After formation, 0.5 to 1.0 mm diameter bubbles 
are injected into the surrounding freestream and reach the freestream 
velocity within a short distance. They are then carried by the air stream 
through the test section. The bubbles are injected at a fairly even rate. 

It is estimated by Hale et al. (73) from photographs that this bubble pro- 
duction rate is as high as 500 bubbles per second. 

The bubble generator console is placed in a convenient location near 
the tunnel. Three lengths of plastic tubing carry the compressed air, 
helium, and bubble film solution from the console to the injector. The 

•> 
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console contains (nlcro'-meterlng valves to control the constituent flow- 
rates, a pressurized cylinder to hold the bubble film solution, and inlet 
connections for high-pressure air and helium. 

An American Optical Co. slide projector, model 1389, served as the 
light source for illumination of the bubbles. A vertical sheet of light 
U cm thick was produced to illuminate the rx-plane of the bubble flow 
pattern. This was provided by using a cardboard slide which is opaque 
except for a thin slit cut out for light passage to the test section. 

Tri-X Pan, a very light-sensitive photographic film rated at ASA 
i«0C, along with a large camera aperture of f 2.0, were employed in order 
to obtain photographs of acceptable contrast. Moreover, the film was 
substantially overdeveloped to compensate for underexposure; this is 
equivalent to using a severely light-sensitive film rated at approximate- 
ly ASA 6,000. 

A typical procedure beoins with adjustment of the flow metering 
valves for air, helium, and bubble film solution in order to obtain a 
high generation rate of small bubbles. Then a Polaroid Land Camera, model 
l80, is used to obtain preliminary photographs since immediate results are 
obtainable. Once the lighting and camera orientation, camera settings, 
and glare are resolved the desired production photographs are taken with 
a 35 mm Minolta single lens reflex camera, model SRT200, with a itS mm 
lens. 

Each bubble leaves an image of its path, or a streak, on the film 
and this streak becomes a record of the air movement within this time 
frame. The record of a number of streaks on a single photograph produces 
a history of the air motion. The intensity of each streak on the film is 
not dependent on shutter speed, but rather on the aperture and bubble 
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I velocity, as well et the sensitivity of the film Itself. The length of 

I pethllnes end (for the slower shutter speeds) the number of pethlines 

I 

ere, however, strongly dependent on shutter speed. A relatively large 
number of streaks provides an overall picture of the fiow pattern. A 
much smaller number though Is useful to see details of ^he air motion 
that are evident only by following individuai streaks for an appreciable 
distance. Production photographs were typically taken at shutter speeds 
ranging from 1/60 to l/k second for each flow condition. 

k. 3.2 Smoke -W I re 

Although smoke particles may not be observed individually, small and 

I 

discrete smoke streaklines can be located within the flowfield so that 

! 

I small scale details may be studied. For example, zero velocity points 

within a large recirculation zone may be identified. 

This is accomplished with the smoke-wire technique, originally devel- 
oped by Raspet and Moore in the early 1950s. It has been improved and ex- 
tended since then to generate sheets of very fine streaklines (7^, 75, 76). 
However, much of the present smoke-wire capability results from the devel- 
opment of Nagib and his coworkers (77, 73). 

The method consists of a fine wire positioned in the flowfield, coat- 
ed with oil, and heated by passing an electrical current through the wire. 
As the wire is coated with the oil, small beads of the oil form on the 
. wire and at each of these beads the smoke filaments originate when the 

wire is heated. Cornell (79) studied the materials used In generating 
"smoke" for flow visualization, and explained that particles created with 
this technique are actually a vapor-condensation aerosol. That is, they 
, arc vary omall liquid particles near 1.0 x 10 ^ mn In diameter rnthcr 

i 
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than products of a combustion procass. HoMaver, In kaaping with tha lit* 
aratura, thasa particlas ara rafarrad to as smoka. 

Tha apparatus, soma of which is shown In FIgura 12, consists of a 
0.28 mm diamatar and 0.6 m long constantan wira, a woodan scaffold sup* 
port frama, a k cc madical syringa, a 0.4$ kg waight, and a solid state 
powar supply, modal TR36*4 h, manufacturad by Elactronic Rasaarch AssocI* 
atas, Inc. Tha support frama holds tha medical syringa oil rasarvoir on 
tha uppar horizontal member, allowing the wire to pass througn the syringe 
needle, the test section, and finally through the lower horizontal member. 
Tha frame is easily relocated so that local details at virtually any 
axial station can be examined. 

Tha wire Is kept taut as the weight Is attached to Its lower end and 
Is free to move as tha wire elongates from thermal expansion. The elec* 
trical leads of the power supply are connected to each end of the constantan 
wire, and conduct about 3«5 amps at 25 volts d.c. The same slide projec* 

m 

tor described for the soap bubble experiment provides the vertical sheet 
of light to illuminate the rx-plane. 

The smoke*wlre technique requires air velocities of less than approx- 
imately 10 m/s; otherwise ‘‘he oil droplets are blown off the wire before 
it is properly coated, and also the smoke diffuses too quickly. Hence 
the smoke experiments were conducted at lower air flow rates than other 
experiments. 

The procedure benins with applying a slight force on the syringe 
plunger, which causes several oil drops to run down the wire, coating it 
with small beads while the wind tunnel is not running. A uniform coating 
of fine oil beads is necessary, as each bead is a source of a distinct 
streakline. Several types of oil were tried, although highly refined 
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kerosene it utilized for ell production photogrepht. Then the wind tur* 
nel motor it turned on. The photographer waits approximately I second 
before manually activating the shutter while the required heat is suppli- 
ed to the oil. A heavy sheet of streaklines then appears, typically for 
a duration of about 0.5 second. Only slight experience with this system 
is required for effective photography. Further details concerning con- 
struction of the apparatus and photography of the nonswirling flowflelds 
are found elsewhere ( 80 ). 

The smoke-wire has primarily been used In considerably less compli- 
cated flow situations such as the flow over airfoils and bluff bodies. 

The present swirling flows exhibit very high levels of turbulent mixing, 
causing the smoke to diffuse ve<'y quickly. Since only very localized re- 
sults are possible under such fim conditions, this experiment was con- 
ducted at several axial stations for each flowfield under investigation. 

A large camera aperture of f 2.0 was used for most production photo- 
graphs. As for soap bubbles, faster shutter speeds such as l/60 second 
provide a more detailed observation whereas slower speeds, for example, 
1/8 second, give a more time-averaged view. Therefore, shutter speeds 
for each flow condition typically include at least l/60, 1/30, 1/15, and 
1/8 second. The same film development procedure used for soap bubble 
photographs is employed here. 

k.3.3 Tufts 

Perhaps one of the simplest ways of visualizing a flowfield is by 
inserting a wire grid with a tuft attached at each intersection point of 
the mesh. These tufts respond to the flow patterns and generally point 
in the local flow direction. Furthermore, their vascillating motion 


33 


indicate* the degree of turbulence in steady mean turbulent flow situa- 
tions. Very large tuft fluctuations, however, can be produced by the 
combined effects of mean flow unsteadiness and turbulence. Such tuft 
screens have been employed to visualize the vortex shedding behind swept 
wings by Bird (81) and Cersten (82). They have also been used by HcMahon 
et al. ( 83 ) to observe the Interaction of an air Jet and a cross wind. 

Observations or photographs from a position normal to the tuft screen 
indicate the flow direction in the plane of the screen. Since the combus- 
tor flowfleld entails three mean velocity components, a tuft screen was 
constructed for each of the rx- and r6-planes. The rx-screen slides into 
the combustor model with onlv slight wall clearance, to minimize the flow 
disturbance. 

The overall dimensions of the rx-screen in the radial and axial direc- 
tions are 29.7 cm and 1.32 m, respectively. Two pieces of extruded alumi- 
num, with 2.5 cm and 0.6 cm as cross-sectional dimensions, are employed 
as axial frame members, spaced to provide a flow passage radial height of 
28.5 cm. Four 1.8 mm diameter steel rods provide radial support and are 
intermittently positioned axially. Music wire of 0.33 mm diameter tra- 
verses the frame with a uniform grid spacing of 1.9 cm radially. A uni- 
form axial spacing of 1.9 cm, 3*2 cm, and A. 4 cm is used for each of 
three axial sections of length 6l cm, i«0 cm, and ^6 cm, respectively. 

The tufts consist of l.l mm diameter cotton thread and are 16 mm in 
length. In construction, the tufts are consistently tied around each 
wire intersection point and a 3.0 mm diameter rod whir- was then removed, 
resulting !n an oversized loop around each wire intersection giving en- 
hanced freedom of movement. 
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The rO'tcrcen Is composed of similer meteriels , except thet It try--* 
is « strip of 24 gauge steel 2.S cm wide, wrapped into a 29.7 cm diameter 
ring. It has a uniform grid spacing of 1.9 cm in both directions. The 
tufts are identical for both screens. 

The light source previously described for soap bubble visualization 
was used as well as the same photographic equipment and technique of over* 
developing the film. Again, slower shutter speeds indicate the tlcre* 
averaged motion, and speeds of 1/2, 1/15, and 1/60 second were typically 
used. 


4.4 Mean Velocity Heasuremsnts 
4.4.1 Instrumentation 

Many instruments are used for separately measuring the magnitude and 
direction of fluid velocity (84). However, there are only a few instru- 
ments capable of simultaneously sensing both magnitude and direction. 

The simplest of these senses pressure at certain locations on the surface 
of a sphere, hemisphere, or some other shape of probe. Spherical pitot 
probes using five holes were used by Janes (85) as early as 1915, Heyer 
and Borren (86) in 1926, and Gutsche (87) in 1931. Pien (88) In 1958, 
further developed the technique upon realizing the fact that the three 
pressure measurements on a great circle of a sphere uniquely determine 
the velocity component in that plane. 

The five-hole pressure probe is recommended by Beer and Chigier (13) 
for mean velocity measurements in turbulent swirling flows and several 
investigators have successfully used it (23, 89, 90). The probe employed 
in this study is model DC-I25-12-CD from United Sensor and Control Corp. 
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It Is shown sch«m«t icsl ly In Figure 1} end hes e 3>2 mm diemeter steel 
sensing tip end sheft contelnlrig five tubes. The sensing heed Is hook* 
sheped to el low prcoe sheft rotetlon without altering the probe tip loca* 
tion. At the opposite end of the sheft, which Is 30. S cm long, Is a 
manifold for the five tubes. 

There are basically three methods of operations. The first and most 
direct Is to adjust the orientation of the probe until both pairs o' 
opposing, dl rect Ion-sen* I ng pressures are nulled, which means the p'-obe 
Is aligned with the local flow direction. With the second method, the 
probe Is fixed and the flow direction Is determined from the calibration 
relationship between probe pressures and flow direction. The th I rd method 
is simply a hybrid of the first two, wherein the yaw angle in the horizon- 
tal plane is nulled and the pitch angle in the vertical plane is correlat- 
ed. 

The first method is avoided since complicated probe orientation 
equipment is required. Further, this equipment would probably disturb 
the flowfield more seriously than the other methods. The second approach 
is completely unacceptable. This is because the probe cannot properly 
sense pitch or yaw angles relative to its axis greater than approximately 
60 degrees, and the yaw angle often exceeds this limit. This limit is 
caused by flow separation at the probe tip. Therefore, the third method 
is exclusively employed in this investigation, as pitch angles are always 
less than 60 degrees except in a very small fraction of the measurement 
locations. 

The instrumentation system, in addition to the five-hole pitot probe, 
consists of a manual ir.' erse mechanism, two five-way ball valves, a very 
sensitive pressure transducer, a power supply, and an integrating volt- 
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meter. The probe Is mounted verticelly on the test chember, as shown in 
Figure ]k, using e manual traverse mechanism, model CIOOO-18 from United 
Sensor and Control Corp. This mechar ism, shown in Figure IS, is made en- 
tirely of steel with a linear vernier accurately readable within t0.2Smm. 
This allows the probe to vertically traverse across the chamber radius 
with the capability of manually rotating the probe about its shaft axis 
to null the yaw angle felt by the pressure sensing tip. This yaw angle 
is read from the rotary vernier of the traverse unit which is accurately 
readable within tO.2 degree. 

The differential pressure transducer is model S900 from Datametrics, 

3 2 

Inc. It has a differential pressure range from 0 to 1.3x10 N/m . The 
sensing element is a high precision stable capacitive potentiometer in 
which the variable element is a thin, highly prestressed metal diaphragm. 
Power for the transducer is provided by a model 72IA Hewlett Packard 
power supply. 

The pressure transducer output is read as the d.c. signal from the 
TSI model 1076 integrating voltmeter. Use of an integrating voltmeter 
removes pressure fluctuations from the vibrating tygon tubing connecting 
the probe to the valves and transducer. The five-way ball valves are 
model SS-43ZF2 by Whitey Valve Co. These valves along with the pressure 
transducer are mounted on a small portable platform, shown in Figure 14, 
in order to minimize the length of tygon tubing for Improved dynamic re- 
sponse of the measurement system. Thus the traverse unit and platform 
are moved together to each successive axial station. 

Finally, auxiliary equipment is used, including model 631 -B strobo- 
tac from General Radio, Inc., to measure fan speed. Also, a barometer/ 
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thermometer unit from Cenco Corp. Is used for local pressure and temper' 
aturci readings. 

k.k.2 Measurement Technique 

The basic irndsurement technique entails aerodynamical ly nulling the 

yaw In the horizontal plane by rotation of the probe about Its vertical 

shaft and then reading two differential pressures (p^ - p^ and p^ - p^) . 

These pressures along with the yaw angle B are used In data reduction to 

obtain pitch angle 6 In the vertical plane and the three mean velocity 

components. The data reduction employs two calibration curves which were 

obtained for a single calibration velocity. The underlying principle is 

that the calibration is independent of probe Reynolds number Re , which 

P 

is based on probe tip diameter. 

Careful calibration experiments reveal that this condition exists 

for Re i 1090 , or a local velocity of 5.^ m/s. Hence measurements of 
P 

such low velocities suffer from a necessary calibration error. However, 
this error affects the velocity measurements typically by less than 6 
percent for Re^ 2 ^00, corresponding to a local velocity greater than 2.0 
m/s. 

Prior to production measurements, the five-hole pitot rotary vernier 
must be zeroed for yaw so that the x and 6 axes of the measurement coordi 
nate frame coincide with those of the test section, which are illustrated 
In Figure lA. This Is accomplished by rotating the probe until the yaw 
is aerodynamical ly nulled near the center of the test section inlet for 
norswirlirg flow. The radial axes of the probe and the test section coin 
clde since the traverse unit base has been carefully machined to mount 
vertically on the plexiglass tube. 
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Each traverse begins at the centerline, and this location is verified 
if the probe casts its shadow in the center of a laser beam which passes 
upstream along the center of the test section. If necessary, the traverse 
mechanism is carefully shimmed until the above condition is obtained. 

After the pressure transducer is zeroed, the measurement procedure 
for each location within a traverse begins with rotation of the probe un- 
til the yaw is aerodynami cal ly nulled. This is indicated by a zero read- 
ing for Py"P£i where the pressures are identified in Figure 13. The 
resulting probe rotation angle <|/ is read from the rotary vernier. Then 
the five-way switching valves are set so that P,, ~ P. is sensed by the 
transducer. Finally, the reading of Pq * P y similarly obtained. 

Hass flow rate entering the test section must be determined for the 
calculation of inlet Reynolds number Re. and spatial mean inlet velocity. 

U 

This quantity is obtained for ea<.h flowfield from a measurement of dyna- 
mic pressure with a conventional pitot-static probe. The probe tip is 
located on the centerline and slightly upstream of the contour nozzle 
exit. It has been experimentally verified with a radial traverse of the 
nozzle discharge that, in this facility, the centerline exit velocity is 
extremely close to the spatial mean velocity. This results from the fact 
that the boundary layer momentum deficit is counterbalanced by a slightly 
h i gher-than-average momentum value just outside the boundary layer. Thus 
the nozzle exit spatial mean velocity for each flowfield is obtained from 
a single pitot-static reading. Finally, a continuity correction for the 
slight flow area difference between the measurement and exit stations pro- 
vides the necessary value. 

The ambient pressure and temperature were monitored near the facility 
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for the determination of density. Also, the wind tunnel fen speed was 
frequently checked with the strobotac. 

An initial experiment was performed to investigate the asymptotic in- 
variance of Re . on the flowfield. This is needed to ensure that the 
d 

dimensionless measurements will be applicable to actual combustor hard- 
ware. which often operates at considerably higher Re^ valuts than those 
desired for the present test fa^'ility. Measurement traverses were per- 
formed at selected axial stations at increased values of Re^ for each 
flowfield. Adjustment of the varidrive fan motor provides excellent and 
stable control of Re^. The resulting velocity profiles are discussed in 
Chapter V I . 

i«.4.3 Data Reduction 

The differential pressure readings from the five-hole pitot probe 
are utilized directly to obtain the square of the vector velocity. In 
turbulent flow conditions, the time-mean pressures give rise to the time- 
mean average of the square of the vector velocity V^. Since 

7^ . + 7^ (4.1) 

where V is the fluctuating portion of the velocity magnitude and tne 

overbars denote time-averaging, It is slightly incorrect to infer that 

is equal to the square of the magnitude of the time-mean velocity vector 
-2 — T 

V . However, the fluctuation term V'^ is not known and very little in- 
formation is available for the effect of turbulen>'.e tn swirl flows on 
pressure probes, which is probably considerable for turbulence intensi- 
ties greater than about 20 percent (12). Furthermore , tha procedures for 
making corrections for turbulence levels are long and tedious, and even 
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then the confidence in their applicability is unknown. Therefore, no 
attempt is presently made to incorporate such corrections, and the deduc- 
ed velocity is taken to be the time'mean velocity magnitude V which is 
written without the overbar from here onward. 

The data at each measurement location are reduced with the FORTRAN 
computer program listed in Appendix C by first calculating the pitch co- 
efficient (Pf^ ■ P 5 )/(Pq ■ Py) • From this value a cubic spline interpola- 
tion technique is used to obtain the pitch angle 6 in the vertical plane 
from the calibration characteristic presented in Figure 16. The result- 
ing value of 6 is similarly utilized to deternine the velocity coeffi- 
cient pV /[2(p^ - p^) ] from the calibration cnaracterist ic given in 
Figure 17. 

Values for V as well as the axial, radial, and swirl velocity com- 
ponents u, V, and w, shown in Figure I8, are easily calculated from the 
velocity coefficient, pi*ch angle, and yaw angle 3, which is in the hori- 
zontal plane. The latter angle is given by 

3 - 360» - K, (i*.2) 

where is the probe rotation angle read on the rotary vernier of the 
traverse mechanism. The magnitude of the velocity vector is given by 


V 





(4.3) 


and the velocity components are obtained from 


u ■ V cos 6 cos 3 
V ■ V s i n6 


and 


(4.4a) 

(4.4b) 
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w ■ V co«A ting 


(i*.4c) 


Further calculations Include that for mass flow rate and axial flux 
of angular momentum, which are respectively expressed as 

R 

m ■ 2n / p ur dr (^.5) 

c 

and 



These are evaluated from the velocity measurements with the data reduc* 
tion computer program. 

Cal ibrat ion 

Proper calibration of the five-hole pressure probe is extremely 
critical to the accuracy of the experimental work. Consequently, consid- 
erable care was exercised in the construction and use of the calibration 
equipment. The calibration equipment consists of a small air jet, a ro- 
tary table, a probe mounting bracket, and the instrumentation system pre- 
viously described. The calibration Jet supply line consists of a com- 
pressed air line, which delivers the desired flow rate through a small 
pressure regulator and a Fischer and Porter model IOAI735A rotameter. 

The jet housing consists of an effective flow management section follow- 
ed by a contoured nozzle with a 3*5 cm diameter throat. 

The rotary table is model BH-9 from Troyke Manufacturing Co. Its 
rotary vernier is readable within ±0.5 minute. As shown in Tigure 19, 
the aluminum probe mounting bracket is secured to the rotary table, and 
it supports the probe which rests in a cylindrical steel collet. 
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The motion of the rotary table orients the probe at the desired pitch 

angle, whereas the yaw is aerodynamical 1y nulled. The probe sensing tip 

remains at the centerline within the potential core of the Jet and less 

than one throat diameter downstream of the nozzle discharge plane. The 

pitch angle is zeroed by very carefully aligning the probe shaft in a 

plane parallel to that of the nozzle throat. 

The calibration procedure consists of recording the voltage output 

from the pressure transducer for differential pressures P|,j " Pj ■nd 

Pc ~ Py* where these pressures are identified in Figure 13. These data 

are obtained at 5 degree increments in 4 over the range -55* s 6 s 55". 

The measurement technique requires the entire calibration to be conducted 

at a constant Jet velocity. This is permitted since the dimensionless 

calibration coefficients are Independent of Re for Re > 1090 (5.^ m/s) 

P P 

determined by careful calibration experiments. Measurements of much 
lower velocities suffer from a calibration error. However, this is rare- 
ly over 6 percent for Re^ > 400 (2.0 m/s). 

Figures 16 ana 17 show the calibration characteristics from which 6 

2 

and velocity coefficient pV /[2(pc-Py)] are obtained, respectively. 

Both curves exhibit considerable symmetry, as the five pressure sensing 
holes are almost symmetrical about the probe tip axis. 
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CHAPTER V 


GROSS FLOWFIELD CHARACTERIZATION 

Recirculation zones are important to combustor designers because most 
of the burning occurs in these regions and they exhibit the highest tempera- 
tures (47). The size and location of these regions is determined by flow 
visualization experiments and predicted by numerical computations. The 
predicted time-mean streamline patterns are compared with corresponding 
artistic impressions deduced from flow visualization photographs of tufts, 
smoke, and bubbles responding to the experimental flowfield patterns. Se- 
lected photographs of the visualization experiments are also presented and 
discussed in the subsequent sections. 

5.1 Artistic Impressions of Streamline Patterns 
From Flow Visualization 

Photographs of each of the six flowfields resulting from g ■ 0, 45, 
and 70 degrees with a ■ 90 and 45 degrees have been examined in detail for 
each of the three flow visualization methods currently employed. The 
characteristics of the overall flowfield are illustrated and discussed via 
the resulting time-mean dividing streamline patterns. These are sketched 
in Figures 20 and 21 from information obtained from the entire collection 
of flow visualization photographs. Results from the smoke-wire experiment 
are utilized near the inlet, whereas tuft and bubble data are used in 
approximating the size and shape of the recirculation zones downstream. 



Bubble flow petterne alto reveal the existence of a precessing vortex 
core, which occurs downstream of the central region. 

The resulting streamlines for the three cases of the sudden expan- 
sion a ■ 90 degree geometry are shown in Figure 20. The nonswirling flow 
sketch in part (a) exhibits a large corner recirculation zone which is in 
excellent agreen«nt with the corresponding streamlines from the measure- 
ments of Chaturvedi (15). Appropriate prior measurements for the present 
swirling flows have not been found. For the moderate swirl vane angle 
case of ^ degrees a central recirculation region appears in conjunc- 

tion with a decrease in size of the corner zone. A thir. precessing vor- 
tex core, discussed by Syred and Beer (91), is observed near the center- 
line extending from the end of the central region to the test section 
exit. The axial location of the upstream varies, ranging approximately 
from x/0 ■ 1.25 to 1.75. This vortex core is essentially a three- 
dimensional time-dependent instability which occurs as a swirling region 
of negligible axial velocity whose center winds around the test section 
centerline. A further increase In vane angle to <(i ■ 70 degrees results 
in continued enlargement of the central zone; however, the corner bubble 
Is essentially unaffected. For this flowfield the vortex core is expand- 
ed in the radial as well the axial direction, as its upstream edge fluc- 
tuates from approximately x/0 ■ 1.0 to 1.5. 

The corresponding sequence of dividing streamlines is found in Fig- 
ure 21 for the gradual flow expansion case with o ■ kS degrees. As with 
Figure 20, the nonswirling flowfield exhibits exce I lent agreement wi th cor- 
responding measurement of Chaturvedi. No evidence of a corner zone is found 
from examination of the photographs for the intermediate vane angle case, 
although the flow pattern is otherwise very similar to the corresponding 


flow for 0*90 degrees. The flowfleld with ^*70 degrees produces exact- 
ly the seme effects es in the abrupt expansion case of Figure 20. 

5.2 Flow Visual Ization Results 


5.2.1 Tufts 

Tuft visualization is very important in that it supplies an overall 
view of zones very clearly. As discussed In section 4.3.3« photographs 
at various shutter speeds were obtained. Slower speeds show more of the 
temporal behavior, although the tufts are sometimes not distinctly visi- 
ble in portions of the flowfield. Some of the more noteworthy photographs 
are presented here, encompassing a range of shutter speeds. Velocities in 
recirculation zones are often somewhat lower than in other portions of the 
flowfield, and thus under such conditions there may be insufficient drag 
on a tuft to align it accurately with the local flow direction. However, 
this is taken into consideration in interpreting the photographed results. 

The three sudden expansion flowfields are characterized in Figure 22 
as photographs of the rx-screen tufts. Figure 22(a) identifies the corner 
zone reattachment point for the nonswirling flowfield. It is found from 
the length scale as y./D ■ 2.1, wherein tufts showing no air motion are 
considered stagnation points. Further, an indication of the turbulence 
level is depicted as the rather slow shutter speed of 1/15 second was 
used. 

For swirl vane angle $ ■ 45 degrees. Figure 22(b) is a photograph 
taken at 1/125 second. The swirl is in the counterclockwise direction 
when viewed from downstream. Both the corner rnd central recirculation 
regions are clearly visible, and their downstream stagnation points appear 
to be at x/D >0.4 and x/0 “ 1.6, respectively. The large vane angle case 
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is photcKir jpi>ed in Figure 22(c), again at I/I2S second shutter speed. Al- 
though the upper corner Is not mcII illuminated in this photograph, the 
lower corner region Is seen to extend to x/D ■ O.A and the central zone to 
x/D - 1.75. 

Photographs of the flowfields with a ■ AS degrees are shown in Figure 
23 . The nonswirling case in Figure 23(a) Is almost Identical to Figure 
22(a), with a corner zone length of x/D ■ 2.0. For ^ ■ AS degrees. Figure 
23 (b) shows no evidence of a corner zone using I/I2S second shutter speed 
and the central region apparently extends to x/D • I.S. Finally, the ^ ■ 

70 degree case exhibits a central region extending downstream as far as 
x/D • l.8S for l/IS second shutter speed. 

Figures 2A and 2S are photographs of the re-screen of tufts for a ■ 

90 and • AS degree flow conditions at x/D - O.S, 1.0, I.S, and 2.S. The 
shutter speed is 1/2 second In order to capture information regarding tur- 
bulence levels as well as flow direction. The camera is positioned slight- 
ly downstream and on the centerline of the test facility, and hence an 
annular ring of distorted tufts appears around the tuft screen because of 
the reflectivity of the plexiglass tube. The variation in brightness of 
the tufts Is merely caused by the angle at which the light strikes each 
particular tuft. 

Figure 2A(a) apparently shows noticeable but random asymmetry in tur- 
bulence levels at x/D - 0.5, whereas Figures 2A(b) and 25(a) seem to re- 
veal a horizontal strip of highly fluctuating tufts. However, it may be 
unrealistic to conclude from this that such a region of intense turbulence 
'exists because the effect of gravity acting on tufts pointing vertically 
upward or downward is different from that on horizontal tufts. Figure 
25 (b) shows that at x/D - 2.3, large direction fluctuations occur only 
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neur th« centtrHne. This probably results from the tendency of the cen* 
ter of the precessinQ vortex core to randomly oscillate about the test 
section centerline (91). 

S.2.2 Smoke -Wire 

Local details in the nonswirling flowfields are clearly revealed 
through the visualization of streaklines indicated from the generation of 
illuminated smoke, as discussed in section ^.3.2. In the swirling flows 
strong mixing diffuses the smoke so that streaklines are not distinguish- 
3 ble. However, under such conditions recirculation zone outlines are 
visible, especially in the region near the smoke- gene rat ion wire. A 
selected photograph is exhibited and discussed for each of the six flow- 
fields considered in previous sections. 

The corner recirculation bubble in the nonswirling case with wall ex- 
pansion angle a ■ 90 degrees is revealed in Figure 26(a) using 1/IS second 
shutter speed. Also, the radial location of the zero velocity point with- 
in the upper and lower corner bubble is estimated to be approximately 
0.15 0 from the respective walls of the test facility. This agrees with 
the velocity measurements presented in Chapter VI. 

The moderate swirl vane angle flow with a 1/30 second exposure is 
shown in Figure 26(b). The shortened corner zone is easily identified 
because the adjacent downstream flow contains no smoke near the inlet. 

This bubble is seen to extend to approximately x/D * O.kS. The tufts in 
Figure 22 indicate a slightly shorter zone ending at approximately x/0 ■ 
O.it. The upstream portion of the central zone is also clearly seen, as 
low velocity fluid carries a dense mass of smoke which slowly moves up- 
stream of the inlet. Further, the precessing vortex core is seen to 
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contain no kmoke. Sine* It exhibits negligible exiel velocity, the voiokc 
Is essentially carried around It by the high velocity fluid outside the 
core. 

Figure 26(c) is a photograph at 1/8 second of the # ■ 70 degree case 
wherein this core is not as distinct. In this case some smoke has diffus* 
ed into the core due to a slightly longer delay before activating the cam- 
era shutter. Observe that both the corner and central zones near the 
inlet reveal that the only change from those for t ■ 4S degrees is a 
si Ight ly wider central region. 

The same swirl sequence for the gradual expansion geometry with a ■ 
degrees Is presented In Figure 27> The nonswirling case at 1/15 second 
reveals about the same details as the corresponding case in Figure 26. 

The 4 • ^5 degree flowfield at 1/60 second exhibits no evidence of a cor- 
ner bubble and a thinner central bubble at the wire position of x/D ■ 

0.25* The existence and location of the vortex core is very similar to 
that in Figure 26(b). The large vane angle case shown in Figure 27(c) 
using 1/15 second reveals no evidence of a corner zone. Also, the central 
zone at x/D • 0.25 appears considerably wider than that for ^ ■ 45 degrees 
shown in Figure 27(b). 

5.2.3 Neutral ly-Buoyant Soap Bubbles 

Soap bubbles injected into the air flow upstream of the test section 
trace pathlines clearly when illuminated. This technique is discussed in 
detail in section 4.3.1. Photographs of illuminated bubbles are taken 
with various shutter speeds and camera positions. In relatively lower 
turbulence intensity portions of the flowfield mean flow directions can 
be obtained by ensemble averaging local tangents to pathlines traced out 
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by so«p bubbles. This h«lp» dafinc the flowfljld geometry In lermt of 
the outline of rec irculet ion regions. 

A sample flow visual izat ion photograph is presented in Figure 28(a) 
corresponding to the zero swirl, 90 degree expansion angle flowfield. 

The photoaraph, taken with a relatively long time exposure (1/8 second), 
clearly shows a great number of individual oathlines. Photographs of 
this type can be used to indicate regions of highly turbulent flow such 
as that near the centerline of the flowfield which exhibits smoother, 
stralghter pathllnes. In addition, the outline of the corner recircuia* 
tion region can be estimated from Figure 28(a) (and numerous additional 
photographs taken at the identical run condition). For this geometry the 
mean stagnation point defining the end of the recirculation zone appears 
to be at x/0 ■ 2.0. 

A photograph with ^ ■ kS degrees and 1/8 second shutter speed is 
shown in Figure 28(b), where the precessing vortex core is clearly seen 
extending fr.om x/D ■ 1.5 to the exit, its upstream edge fluctuates ran- 
domly from approximately x/0 ■ 1.25 to 1.75. The corner bubble is observ- 
ed in both the upper and lower portions of the flowfield, extending to 
approximately x/0 ■ 0.i«, which agrees almost exactly with the smoke flow 
pattern for this flowfield in Figure 26(b). A photograph using 1/15 sec- 
ond shutter speed is presented in Figure 28(c) for the ^ 70 degree flow- 

field where a thicker vortex core is obvious. Also, It extends farther 
upstream to x/0 * 1.25 in this photograph, although this fluctuates from 
about x/0 " 1.0 to 1.5. The corner zone is faintly visible here, and its 
axial length seems to extend to about x/0 ■ O.k. 

Figure 29(a) is a 1/15 second exposure of the nonswirling case with 
a ■ A5 degrees. Essentially the same result is found as in the 
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corresponding a ■ 90 degree flowfield in Figure 28(a). The intermediate 
vane angle case is visualized using 1/^4 second shutter speed in Figure 
29(b). As with the smoke flow pattern in Figure 27(b), there is apparent- 
ly no evidence of a corner bubble. Otherwise, there is little difference 
uetween Figures 29(b) and 28(b). The ♦ • 70 degree flow shown In Figure 
29(c) is photographed at 1/8 second, and agrees with the smoke pattern in 
Figure 27(c) that there is no apparent corner zone. Further, the vortex 
core extends upstream to x/D - 1.25, as in the corresponding swirl case 
in Figure 28(c) . 


S.3 Predicted Streamline Patterns 

The inlet boundary condition for the flowfield prediction program 
consists of a flat velocity profile as in plug flow for both u and w, 
whereas v ■ 0. This approximation has often been employed when pertinent 
measurements are not available (46, 47). It was discovered near the end 
of this investigation from the five-hole pitot velocity measurements for 
the swirling f lowf ields, that in the test facility, the Inlet velocity pro- 
files are actually much different In shape and magnitude from those assum- 
ed for the general predictions. These measurements are presented and fur- 
ther discussed in Chapter VI. This difference results from: (a) Imperfect 

blade efficiency, (b) the existence of a hub, and (c) the fact that the 
downstream edge of the swirl vanes of the test facility is actually locat- 
ed approximately 4 cm upstream of the flow expansion corner, where x/D “ 

0. This swirler location allows the central recirculation zone to begin 
upstream of x/0 » 0, thereby changing the velocity profiles there. Hence 
swirling flow streamline predictions in Figures 30 and 31, parts (b) and 
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(c) , «r« not expected to exhibit excellent egreement with experiments, 
perticulariy near x/D ■ 0. 

Figures 30 and 31 are calcuiated and piotted by computer to show time* 
mean streamline patterns for a ■ 90 and 1*5 degrees, respectively. Although 
the solution of the swirling flows is not entirely realistic, the trends 
are generally correct. The nonswirl predictions do not suffer from the 
inlet boundary condition error described above. The results shown in Fig- 
ure 30(a) exhibit excellent agreement with the flow visualization results 
in Figure 20(a) and also with the streamline pattern resulting from corre- 
sponding measurements of Chaturvedi (15). The intermediate case of ^ ■ 
k5 degree vanes shows reasonable agreement with the results in Figure 
20(b) except at x/0 ■ 0. Discrepancies include central recirculation zone 
lengths of 1.35 0 and 1.6 D for predictions and visualization, respective- 
ly. The corner recirculation rei^l'Vi lengths, however, are in exact agree- 
ment. Surprisingly, for ^ ■ 70 degrees the predicted and photographed 
central zone lengths are x/D ■ 3.0 and 1.75, respectively. Further, no 
corner zone is predicted, whereas the photographs reveal one with axial 
length x/D ■ 0.4. 

As with the a • 90 degree case, the nonswirling flowfield with a ■ 

45 degrees, shown in Figure 31(a), is very well predicted by comparison 
witii Chaturvedi 's measurements and the present photographs. No corner 
zone is predicted in this geometry for either the intermediate or large 
vane angle flows, which agrees with visualization results of Figure 21(b) 
and 21(c). For $ ■ 45 degrees, the predicted central zone length Is x/D 
- 1.35wherea that photographed is approximately 1.5. The central region 
length for the strongly swirling flow from computations and visualization 
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is x/D ■ I.7S and I.SS. respectively. Hence the central zone is more 
accurately predicted for a ■ 45 degrees than for a ■ 90 degrees, 

5.4 Parametric Effects 

Streamline plots, including those deduced from flow visual izat ion as 
well as those predicted, allow recirculation zones to be characterized 
parametrically. The parametric effects are illustrated In Figures 32 and 
33 for predictions and visualization, respectively. Figure 32 shows the 
predicted effects of a and ^ on the corner and central recirculation zone 
lengths. Several observations should be noted. First, zone lengths are 
only slightly affected by a, as found previously for the corner bubble 
under isothermal nonswirling conditions (15, 31). Further, at these high 
Reynolds numbers, the corner region length is independent of Reynolds num- 
ber (33. 34). As a decreases from 90 to 45 degrees, the zone lengths tend 
to decrease slightly and the inlet flow is encouraged to impinge on the 
confining walls. Second, the corner bubble length decreases upon increas- 
ing from 0 to 45 degrees, and also the central region and the precessing 
vortex core appear. Increasing the vane angle to 70 degrees slightly en- 
larges the cen^'al zone and the vortex core. 

Visualization results exhibit close agreement with predictions for 
the shape of be ‘ h corner and central regions for ^ - 0 and ^ ■ 45 degrees. 
Further, this also occurs for ^ > 70 degrees with the gradual expansion 
geometry where a ■ 45 degrees. However, for the abrupt expansion case 
where a - 90 and ■ 70 degrees, the visualized corner zone is not predict- 
ed. Moreover, the central region length is seriously overpredicted at x/D 
■ 3.0 as compared with 1.79. A lack of prediction accuracy for swirling 


53 


flows Is expected because of compliceted inlet conditions and/or tbe lack 
of realism in the turbulence model; however, it is surprising that the a ■ 
4S degree case exhibits close agreement while that for a ■ 90 degrees does 
not . 


I 
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CHAPTER VI 


VELOCITY MEASUREMENTS, PREDICTIONS, AND 
COMPARISON WITH EXPERIMENT 

Each of the six flowfields under investigation is first character- 
ized in terms of axiai and swiri time-mean velocity measurements obtained 
using a five-hole pitot probe. Then generai velocity predictions using 
estimated inlet conditions are presented. A special test case is do- 
scribed also which serves as a validity check of the computer code for 
the abrupt expansion angle a ■ 90 degrees using measured inlet conditions. 
Next, reliability of predictions and measurements is discussed, and final- 
ly, the results of another validity test for a non-expanding geometry is 
described. 


6.1 Five-Hole Pitot Probe Velocity Measurements 

Details of the five-hole pitot measurement apparatus and technique 
are discussed in section k.k. Readings were taken at approximately 20 
positions across the test facility radius at the various axial stations 
for which velocities are presented. Before production measurements be- 
gan, readings were taken for ^ ■ 0, ^5, and 70 degrees at various values 
of inlet Reynolds number. Comparison of these data determined the 
Reynolds number required to ensure that the nondimensional ized measure- 
ments are independent of Reynolds number. A sample of these data is pre- 
sented and discussed in section G.U.2. 
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Figures 3Si and 36 show the axial and swirl velocity profiles 

for 9 ■ 0, ^5. and 70 degrees with a ■ 90 degrees. Note that the scale 

of dimensionless u/u^ shown at the inlet in Figures 3^(a) and 35(a) is 

different from that of w/u which is found in Figure 35(b). These scales 

o 

apply to all figures through section 6.3. Radial velocities are consis* 
tently much smaller than the axial and swirl components, and they are 
indicated in the dividing streamline patterns deduced from flow visual- 
ization experiments depicted in Figures 20 and 21. Therefore, they are 
not presented here. Figure 3^ shows the corner recirculation region for 
the nonswirling flow on the verge of reattachment at x/D ■ 1.5 and also 
at 2.0, whereas it occurs at x/D ■■ 2.1 for the corresponding visualiza- 
tion results shown in Figure 20(a). The corresponding measurements of 
Chaturvedi (15) yielded a value of x/D > 2.3. 

Both of the inlet profiles in Figures 35 and 36 differ in magnitude 
and shape from the fiat profiles used as inlet boundary conditions for 
the general predictions. This Is discussed in section 5.3 where it is 
mentioned that the swirl vanes in the facility are not 100 percent effi- 
cient. It can be seen in Figure 35 that a 45 degree swirl vane angle 
produces a maximum swirl flow angle [tar. ^ (w/u) ] of 30 degrees at the 
test chamber inlet. The maximum swirl flow angle for the 70 degree vane 
angle case shown in Figure 36 is only 34 degrees. Thus there is only a 
slight increase of swirl flow angle, although these two flows are differ- 
ent in that the inlet profiles are considerably more sharply peaked for 
the ^ • 70 degree case. Figures 35 and 36 show zero u and w velocity 
near the axis at the inlet, but actually the probe was insensitive to 
the very low velocities there. This is consistent with the flow 
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visualization result discussed in Chapter V that the central bubble ex- 
tends upstream of the inlet. 

The measurements shown in Figure 35(a) provide no evidence regard- 
ing the existence of a corner zone. This is expected because flow visu- 
alization reveals that the region only extends to x/D > 0.'«. However, 
there is clear evidence of a rather large central zone whose length is 
similar to that shown in Figure 20(b). Although the axial velocity pro- 
files are beginning to flatten in the downstream direction, they retain 
a zero velocity value on the axis, which is consistent with soap bubble 
flow patterns as seen In Figure 20(b). The early erratic behavior shown 
by the swirl velocity profile in Figure 35(b) quickly transforms, exhib- 
iting a solid-body-rotation core with a rather flat profile outside this 
region. 

Figure 36 also reveals the discrepancy regarding inlet velocity pro- 
files. The large central recirculation region causes the downstream 
flow to be accelerated near the top wall. The early erratic behavior 
found in the swirl velocity profile at x/D - 0.5 quickly develops into a 
shape similar to that seen in Figure 35(b). The processing vortex core 
.TOtion discussed in Chapter V results in poor measurement repeatability 
which promotes the irregular behavior within this core region. 

Figures 37, 38, and 39 exhibit velocities for the same sequence of 
flowfields with a ■ 45 degrees. The inlet profiles were not measured in 
this geometry because the presence of the expansion block interferes. 
Effects of (|i on velocities, similar to those found for the sudden expan- 
sion cases, are found in the flowfield sequence for this test section 
geometry. The niajor difference is that the sloping wall encourages the 
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Inlet floM to accelerate near the top wall. Also, It tends to shorten 
or obliterate the corner rec I rculat Ion region. 

6.2 General Velocity Predictions 

All results are obtained via a nonuniform grid system due to the re- 
sulting enhancement of solution accuracy. Nineteen cells are employed 
In the r-directlon, and they are clustered near the shear layer region 
and along the wall and centerline. The cells are gradually expanding In 
the streamwise direction, and from 21 to 33 are employed as required to 
produce the desired a. The Inlet profiles of axial velocity u and swirl 
velocity w are Idealized as "flat" (I.e., constant-valued) with 100 per- 
cent efficient swirl vanes. Also, It Is assumed that v ■ 0. The exit 
Is always located well downstream of the central zone where the flow Is 
considered parallel with zero radial velocity. Further details are 
briefly discussed In Chapter III. For comparison with other results. It 
should be noted that $ and swirl number S are related by S ■ 2/3 tan 
so that vane angles 43, 60, and 70 degrees, for example, correspond to S 
values of 0.67, 1.15, and 1.83, respectively (12). 

Predicted time-mean axial and swirl velocity profiles for the slde- 

wa’I expansion angle a ■ 90 degrees are shown in Figures 40, 41, and 42 

for swirl vane angle $ ■ 0, 45, and 70 degrees, respectively. The non- 

swirling case In Figure 40 exhibits close agreement with measurements of 

Chaturvedi (13) as shown elsewhere (3). Note the scale of nondimenslonal 

u/u shown at the inlet profile; this applies to all subsequent general 
o 

predictions. Further, the centerline velocity exhibits little change in 
the streamwise direction. 
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Figure 41(a) reve4.!» that with « ■ 4S degrees, the mean axial veio* 
city profiles are dramatically changed. A maximum velocity value occur., 
in an annular fashion near r/D ■ 0.35i although a more flattened shape 
quickly develops before x/D ■ 1.0. It should also be noted that the 
boundary layer on th<i outer sidewall Is relatively thin and Is not easily 
seen on the figures using the chosen scale. Note that the scale for w/u 

o 

shown at the inlet in Figure 41(b) Is different from that of u/u . Both 

o 

of these scales are used in all subsequent figures through section 6.3. 
All of the swirl velocity profiles exhibit sol i d-body-rotat Ion behavior 
near the centerline, even near the inlet where a flat profile is a speci* 
fled inlet condition. The radial location of the station maximum for w 
tends to increase with x/D in Figure 41(b). Swirl as well as axial velo- 
city profiles appear to approach those corresponding to swirling flow In 
a pipe as x Increases. 

Figure 42(a) reveals that the strong swirl case provides a much 
larger central recirculation region at x/D ■ O.St which is caused by 
strong centrifugal effects. This promotes a very large foward velocity 
near the wall and v.most eliminates the corner region, although experi- 
ments show it is unaffected. The radial extent of the central recircula- 
tion zone and the velocity near the wall quickly diminish downstream as 
swirl strength is dissipated. The swirl velocity profile in Figure 42(b) 
are similar to those of Figure 41(b), except much stronger magnitudes re- 
sult from the increased vane angle. 

The same sequence for the a ■ 45 degre.** geometry is found in Figures 
43, 44, and 45. The inlet profiles are not shown here since they are the 
same as those of the corresponding cases in Figures 40, 4l, and 42. The 
nonswirling flowfield shown in Figure 43 is essentially identical to the 
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nonswlrling case with o ■ 90 degrees. For ♦ • l«5 degrees. Figure 
shows little effect of a upon comparison with Figure 41. The central 
region Indicated In Figure 44(a) Is slightly thinner at x/0 ■ 0.5. The 
strong swirl flowfield illustrated in Figure 45(a) exhibits a smaller 
backflow region than for a ■ 90 degrees in Figure 42(a), although it dis~ 
plays larger upstream velocities. Also, the peak in the profile at x/D 
■ 0.5 occurs in an annular fashion, whereas this is not seen in Figure 
42(a). Swirl velocities in Figure 45(b) exhibit a slightly larger magni- 
tude near the top wall than that in Figure 42(b), but are otherwise very 
similar. 


6.3 The Specific Test Case 

Due to the discrepancy regarding axial and swirl inlot velocity pro- 
files, the importance of which was discovered near the end of this inves- 
tigati on, another prediction of the a ■ 90 and $ ■ 45 degree flowfield 
is presented. The measured velocities are specified as the inlet boun- 
dary condition, and the results are seen in Figure 46 along with the mea- 
surements. 

This prediction shows that the large peak of inlet axial velocity 
for x/D > 0.5 is located at r/D ■ 0.25 rather than 0.40 as measured. 

Also, a much smaller central zone is predicted than for the general pre- 
dictions. However, no effect of this discrepancy is seen beyond x/D ■ 
1.5. The measured inlet swirl velocity peak broadens very quickly by 
x/0 - 0.5. Also, a solid-body-rotation core is exhibited there, which 
is somewhat similar to the previous predictions shown in Figure 41. Com- 
parison with measurements shows considerably larger velocity magnitudes 
at x/D - 0.5, but lower values downstream. 
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This dIscrapancY naar tha inlat might rasuit from backflow at the 
iniat boundary which could not ba maasurad bacausa of transiant effects 
or tha lack of probe sensitivity for smail velocities. 

6.4 R liability of Predictions and Measurements 


6.4. i Predictions 


The finite difference mash size employed for the flowfield predic* 
tions resulted from a compromise between the irreconcilable alternatives 
of computing cost and flowfield resolution. Several grid systems were 
utilized to determine the effect of grid density on the solution of the 
present flowfieids. The chosen grid is discussed in Chapter III and an 
example is shown in Figure S>. An increase of 60 percent in the number 
of computational cells produces a change of only 0.3 and l.l percent for 
top w.ill and centerline dimensionless axial velocities, respectively. 

Another important aspect of solution accuracy concerns the converg- 
ence criterion. Rather than use the well established fractional change 
criterion, which is a necessary but not sufficient condition for converg* 
ence, a residual source criterion is employed. This residual is defined 
by equations like 


■ ‘U, - 1 - s; 


( 6 . 1 ) 


which has been obtained from Equation (3.1) and which measures the depar- 
ture from exactness (except for the use of current coefficients and 
source) for the variable ^ at the point P. When each of these quantities 
becomes smaller than a certain fraction of a reference value, the finite 
difference equations are considered to be solved. 
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Although inaccur«ci«t resulting from mesh resolution and degree of 
convergence are insignificant, the predictions suffer from inherent in* 
adequacies in the turbulence model. This is tolerated presently because 
these models for complicated swirling flows are not yet available (A9) . 
The k-c model developed for nonswirling flows has been extended to in- 
clude swirl shear stresses; however, an isotropic viscosity remains 
while swirling flows in combustor geo.'ietrics can exhibit anisotropic be- 
havior (26). However, the test case coniisting of the measurements of 
Vu and Gouldin (29) exhibits very encouraging results as discussed in 
section 6.S. 

6.k.2 Measurements 

It is difficult to estimate the accuracy of the fkw visualization 
results. The soap bubbles probably trace the flow very accurately be- 
cause they are neutral ly-buoyant and approximately only 0.5 to 1.0 mm in 
diameter. The tufts observed in the rx-plane of the test section do not 
respond well to the local flow direction in a low velocity region due to 
the small drag forces acting on them. They are con< ’dered very reliable 
in higher velocity regions, however. The smoke particles have much 
greater density than air and consequently they are affected by gravity 
in low velocity areas and by centrifugal forces in strongly swirling re- 
gions. 

Little information is available concerning the effect of turbulence 
on pressure probes in swirling flows (12). However, it is presumed that 
the five-hole pitot probe is accurate within approximately 5 percent for 
most of the measurements. This value may increase to 10 percent as the 
velocity magnitude falls below approximately 2.0 m/s because the probe 


it not v«ry tcntitiv* to tucH Iom dynamic pretsuret. In addition, tha 
calibration it no longar Indapendent of proba Raynoldt numbar. Tha prat* 
tura trantducar It attimatad to hava a maaturament uncartainty of approx* 
imataly 0.2$ parcant for a typical valocity magnitude of S m/t and 0.75 
percent at 2 m/t. 

An axparlmant to determine the extent to which the five-hole probe 

ditturbt tha flow wat conducted. A dummy probe of approximately the tame 

tize and thape wat interted through the lower wall of the tett taction, 

exactly oppotite to tha five-hole probe entering through the upper wall. 

Axial and twirl velocity prof! let with and without the dummy piobe, which 

it t reverted exactly like the meaturing probe, are thown In Figuret A/ 

and kS for x/D - 1.0 and 2.S, retpect i vely , with a ■ 90 and ^ ■ kS de- 

greet. Clearly the pretence of the dummy probe hat virtually no effect. 

Further, tince the dummy probe hat a negligible influence, these figuret 

show evidence of good measurement repeatability. 

Before beginning velocity surveys, an initial experiment was per* 

formed to investigate the asymptotic invariance of Reynolds number Re.. 

a 

This was needed to ensure that dimensionless measurements will be appli- 
cable to combustors which operate at higher Re. values than those desired 

U 

for the test facility. Measurement traverses were repeated at the same 
conditions for Increasing Re^ for each of the flowfields at x/D ■ 2.5. 

This axial station was chosen because a probe location slightly down- 

stream of a recirculation zone was considered most sensitive to Re.; and 

0 

the central rather than corner zone was selected since most of the com- 
bustion and flame stability occur there. A sample of the resulting 
axial and swirl velocity profiles Is found in Figures ^9(a) and 49(h), 
respectively. Note that there Is little change for 68,000 £ Re . s 92 ,000. 

Q 
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Hence ell measurements of the moderate ^ • kS degree floMflelds are par* 

formed at Re, ■ 78,000. Similarly, all data for a ■ 0 and g ■ 70 degrees 
0 

are obtained at Re . ■ I0S,000 and 78,000, respectively. 

0 

Figure SO shows axial and swirl velocity profiles at x/D • 1.0 with 
■ 90 and g ■ 45 degrees. The measurements above and below the axis 
were obtained from probe entry through the upper and lower wall, respec* 
tively. Although the upper and lower profiles reveal a similar shape, 
they exhibit different magnitudes. This shows clear evidence of a slight 
asymmetry in the flowfield. 

Averaging the mass flow rate calculated from each radial traverse 
In Figure 50 results in a value within 2 percent of that measured at the 
test section inlet. Calculation of mass flow rate from each traverse 
Individually yields a discrepancy of approximately 14 percent of the mea- 
sured value. This asymmetry is Indicative of the lack of agreement which 
may be expected if mass flow rate Is calculated from the upper traverse 
only. 

6.5 Prediction of Coannular Swirling Pipe Flow 
6.5-1 The Flowfields Considered 

Due to the rather poor agreement resulting from tiie special test 
case near the inlet of the test section, another validity test of the 
prediction code was undertaken. This case, experimentally investigated 
by Vu and Gouldin (29), has no expansion and the inlet velocity profiles 
do not exhibit backflow. Therefore, it Is considered less difficult to 
predi ct . 

This test consists of both counterswirling and coswirling flows 
from a coannular jet discharging into a straight pipe. A schematic 
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diagram of the flow assembly Is Illustrated In Figure 5I« The assembly 
consists of two major segments: an Inner flow passage and an outer flow 

channel. The Inner passage consists of a circular tube with inner and 
outer diameters of 3>72 cm and 3.8 cm, respectively. A swirl vane assem- 
bly containing 12 equally spaced vanes with a swirl vane angle of 68.5 
degrees is located 11.9 cm upstream of the discharge end of the tube. 

The outer flow channel Is a radial inflow passage with 2k adjustable 
swirl vanes near its inlet. The outer flow turns from the radial to the 
axial direction whereupon It exits in an annular fashion into a IA.6 cm 
diameter pipe where the mixing of the two streams occurs. 

6.5.2 Counterswirling Case 

The inner and outer flows have swirl numbers of O.kS and -0.51, re- 
spectively. A nonuniform fine mesh of size 30 x 24 (Nl ■ 30, NJ ■ 24) 
is arranged to cover a pipe region of length 0.17 m and radius 0.0727 m. 
Inlet profiles of axial and swirl velocities are taken from the measured 

vilues at the location x * 0.002 m (x/D ■ 0.01375), as indicated in Fig- 

3/2 

ure 52. Inlet k and e (■ /l) are specified via 

k - 0.03 u^ (6.2a) 

I - 0.005 D/2 (6.2b) 

where D ■ 0.1454 m. The exit flow is parallel, with zero axial gradients 
for all variables and v ■ 0. 

After a preliminary calculation of 22 iterations and a further 189 
iterations, final convergence is obtained. The criterion is that the 
maximum normalized residuai source sum for each of the variables mass, u, 
V, w, and p is less than 0.004. Radial profile computations of axial 
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and swirl velocities u and w are shown in Figure 33 at the axiai station 
x/D • 0.137 (x ■ 0.02 m) , and in Figure 5^ at the station x/D ■ 0.27^ 

(x * 0.0l« m) . These figures show veiocities normaiized with respect to 
a fixed quantity u^, taken as 2i.6 m/s as an estimate of the outer fiow 
spatiai average axial velocity, and the experimentai data, quoLt.d by 
Srinivasan and Mongia (1*9), have been appropriately transferred to these 
figures for comparison. The two axial stations chosen lie in a central 
zone and show dramatically the differences. 

Notice that, whereas the agreement between swirl velocity profiles 
is very satisfactory at both locations, the agreement between axial velo- 
city profiles is not. in the present work the central zone is predicted, 
but it is not as Idrge or with such large negative veiocities as found in 
the experiment (29). Finally, Figure 55 exhibits the centerline axial 
velocity development as a function of downstream location. The velocity 
predictions at the nearest gridline to the axis are taken as representa- 
tive of axis values (r/D ■ 0) since 3u/3r is zero there. The central re- 
circulation zone is rather well predicted in the present study, although 
the axial velocity magnitudes are all underpredicted. 

6.5«3 Coswirling Case 

The flowfield for this experiment concerns inner and outer flow 
swirl numbers of 0.58 and 0.5^, respectively. Inlet profiles of axial u 
and swirl w velocities taken from the measured values at the location x 
■ 0.002 m (x/D ■ 0.01375), as indicated in Figure 56. Other boundary 
conditions are as taken previously. Radial profile computations of axial 
and swirl velocities are shown in Figure 57 at the axial station x/D ■ 
0.0688 (x ■ 0.01 m) , and in Figure 58 at the axial station x/D ■ 0.2063 
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(x > 0.03 m) . These may be compared directly with experimental evidence 
given in the corresponding figures. Velocities are normalized as in the 
previous counterswiri care. The present predictions are very encourag- 
ing when compared to the experimental data. Finally, Figure 59 illus- 
trates the centerline axial velocu, development as a function of down- 
stream distance. Again, velocity predictions are taken from the nearest 
gridline to the axis. Notice that predictions of the present study com- 
pare quite well with the experimental data, especially for x/D s 0.5. 
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CHAPTER VI I 


CLOSURE 

7.1 Conclusions 

The present research is concerned with a specific sub-problem of 
the complete three-dimensional problem and is restricted to axisymmetric 
geometries under low speed and nonreacting conditions. This is an area 
in which there is a need for more fundamentfil research, particularly re- 
garding recirculation zones where most of tie burning occurs in reacting 
flows . 

Many factors affect the existence, size, and shape of both the cor- 
ner recirculation zone and the central toroidal recirculation zone. A 
major outcome of the current study is the experimental and computational 
characterization of these for six basic flowfields. Parameter varlationi 
which define these flow conditions are side-wall angle a • 90 and kS de- 
grees with swirl vane angle 0*0, kS, and 70 degrees. The size and 
shape of the recirculation bubbles for each flowfield is illustrated as 
an artistic impression deduced from a collection of flow visualization 
photographs of tufts, smoke, or neutral ly-buoyant soap bubbles responding 
to the flow. 

Photographed results shov. that increasing 0 from 0 to 45 degrees 
produces a shortened corner region and the appearance of a central bubble 
typically extending downstream to approximately x/D ■ 1.7. Also, a pro- 
cessing vortex core is observed which scratches from the end of the 
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central region to the test section exit. A further Increase In ^ to 70 
degrees enlarges the central zone and vortex core with negligible effect 
on the corner region in those flowfields where it occurs. 

The effect of n on the nonswirling flows is negligibie; however, a 
decrease from 90 to ^5 degrees apparently eliminates the corner bubbie 
in each swirling flow. This decrease in a also causes the inlet flow to 
impinge more severely on the top wall, where larger axial velocities 
occur. 

A more detailed experiment consists of the measurement of time-mean 
velocity components In the axial, radial, and azimuthal directions using 
a five-hole pitot probe. These measurements generally agree with the 
flow visualization results and provide a more complete understanding of 
eacn flowfield. Also, they constitute a data base allowing validity 
*’esting of flowfield predictions. Further, inlet velocity measurements 
reveal that the estimated inlet velocity profiles employed in the present 
swirling flowfield predictions are unrealistic. This arises from blade 
inefficiency, the presence of a hub, and the fact that the swirl vane 
exit station in the test facility is actually located k cm upstream of 
the test section inlet. This allows the central recirculation zone to 
begin upstream of the inlet where x/D ■ 0. 

The flowfield prediction computer program is employed in a paramet- 
ric study to computationally determine trends of interest to combustor 
designers. The predicted effects of a and on swirling flows are almost 
the same as those determined experimentally, even though the inlet velo- 
city estimates employed consist of flat profiles with 100 percent effi- 
cient swirl vanes located at the inlet with no recirculation. Further 
predictions using measured inlet velocities as boundary conditions were 
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obtained a$ a special test case to allow appropriate comparison with mea- 
surements for the a • 90 and ♦ ■ ^5 degroe flowfield. Poor agreement re- 
sults at x/D ■ O.Sf although it is excellent for x/D i I.S« This dis- 
crepancy may be due to the lack of realism inherent in the k-c turbulence 
model (extended to include swirl shear stresses), inasmuch as the empiri- 
cal constants pertain to nonswirling flows. Another possibility may be 
an inaccurate specification of velocity and turbulence distributions at 
the inlet to the test section. 

7.2 Recommendations for Further Work 

Fundamental research should be continued in several areas. First, 
time-mean velocity measurements should be obtained ai. additional axial 
stations, especial iy in the upstream region, to provide more details and 
allow more extensive prediction validity tests. Second, movies of the 
flow visualization experiments should be obtained and analyzed to charac- 
terize the flow instabilities and temporal behavior of the recirculation 
zones and the trailing vortex core. Third, further details of the flow- 

I 

fields are needed to allow a more complete fundamental understanding. 

This should consist of turbulence quantity measurements, including all 
of the Reynolds stresses. This would allow the deduction of more sophis- 
ticated turbulence models relating turbulent shear stresses with time- | 

mean velocity gradients. 

. i 

I 

I 
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^tgure 2. The Flowfleld Being Investigated 
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Figure The Three Control Volurnes Associoted With 
Points of the Three Grids 
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Figure 6. General View of the Test Facility 
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Figure 7. View of the Swirl Vane Assembly 
From Upstream 



Figure 8. View of the Swirl Vane Assembly 
From Downstream With the In- 
strument for Measuring the 
Vane Angie 
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Figure 9 
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Figure 10. Apparatus For Neutral ly-Buoyant Bubble Flow Visual izatlon Experiment 
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Portable 



Figure I A. Apparatus for Nean Velocity Heasurements Using 
a Five-Hole Pitot Probe 
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Pitch Angle Calioratlon Characteristic for Five-Hole Pitot Probe 







Figure 19. Calibration Apparatus With Five-Hole Pitot Probe 



Figure 20. 
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Figure 21 . Artistic Impressions of Dividing Stream 
lines With Wall Expansion Angle a « 
45* for Swirl Vane Angles: (a) ♦ » 

0 “, (b) « - 45*. and (c) ^ - 70 ” 
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(b) ^ - 45' 



(c) ^ - 70' 

Figure 22. Flow Visualization Photographs of Tufts 
in the rx-Plane With Wall Expansion 
Angle a ■ 90° for Swirl Vane Angles: 
(a) 4) - 0®, (b) 4) - 45', and (c) 4> - 
70' 
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(a) ♦ - 0* 



(b) ♦ - i*5‘ 



Figure 23. Flow Visualization Ph-jtographs c f Tufts 
in the rx-Plane With Wall Expansion 
Angle a ■ for Swirl Vane Angles: 
(a) - 0®, (b) iji - i*5“, and (c) <|i • 

70" 
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(b) x/D - 1.0 

Figure 2k, Flow Visualization Photographs of Smoke -W I re 
Streaklines With V'.'M Expansion Angle a ■ 
90* for Swirl Vane Angles: (a) x/D “ 0.5 

and (b) x/D ■ 1.0 
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(a) x/D - 1.5 



^ = 90 
0 .? 45 
X/D‘: 2.5 
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(b) x/D - 2.5 

figure 25. Flow Visualization Photographs of Tufts 
in the r6~Plane With Wall Expansion 
Angle a ■ 90* and Swirl Vane Angle (> 

■ 45* for Axial Stations: (a) x/0 ■ 

1.5 and (b) x/D - 2.5 
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(•) ♦ - 0 * 



I A 



(c) « - 70" 

Figure 26. Flow Visualization Photographs of 
Smoke-Wire Streaklines With Wall 
Expansion Angie o ■90* for Swirl 
Vane Angles: (a) ♦ - 0* , (b) ^ 

- i«5*, and (c) ^ - 70* 
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U) ♦ - 0* 



(c) ♦ - 70“ 

Figure 27. Flow Visualization Photographs of Smoke* 

Wire Streaklines With Wall Expansion 
Angle a - 45* for Swirl Vane Angles: 

(a) ^ - 0“, (b) ^ • 45*, and (c) * - 

70* i 
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(a) ♦ - 0* 



(b) « - i*S‘ 



(c) 4> - 70" 

Figure 28. Flow Visualization Photographs of Pathlines 
Indicated by Illuminated Neutral ly-Buoyant 
Soap Bubbles for Wall Expansion Angle a ■ 
90“ and Swirl Vane Angles: (a) ^ ■ 0“, 

(b) (> - 45®, and (c) (p - 70“ 
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(b) 4 - i*5’ 



Figure 29. 


(c) ♦ - 70* 


Fiow Visualization Photographs of Path- 
lines Indicated by Illuminated 
Neutral ly-Buoyant Soap Bubbles for 
Wall Expansion Angle a - 45“ and 
Swirl Vane Angles: (a) 4 ■ 0“ , (b) 

4 ■ 45“ , and (c) * ■ 70“ 
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(ii) Corner Recirculation Zone 

f i jUre }2. Predicted Effect of Wall Expan- 
sion Angle a and Swirl Vane 
Angle <(> on: (a) Corner Re- 

circulation Zone Length and 
(b) Central Toroidal Recircu- 
lation Zone Length 
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(b) Central Toroidal Recirculation Zone 
Figure 32. (Continued) 
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(a) Corner Recirculation Zone 

Figure 33. Experimental Effect of Wall Ex- 
pansion Angle a and Swirl 
Vane Angle ♦ on; (a) Corner 
Recirculation Zone Length and 
(b) Central Toroidal Recircu- 
lation Zone Length 
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;ral Toroidal Recirculation Zone 
Figure 33- (Continued) 
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Figure 35. Measured Velocity Profiles for Wall Expansion Angle 
and Swirl Vane Angle 9 » 45“ 








Figure 37. Measured Velocity Profiles for Wall Expansion Angle 
and Swi rl Vane Angle ^ = 0“ 














Figure 4l. Predicted Velocity Profiles for Wall Expansion Angle 
and Swirl Vane Angle ^ ■ 1*5® 
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Figure Uk. Predicted Velocity Profiles for Wall Expansion Angle a • kS 
and Swirl Vane Angle ^ > kS* 
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Figure k5. Predicted Velocity Profiles for Wall Expansion Angle 
and Swirl Vane Angle | > 70* 



Figure 46. Predicted and Measured Velocity Profiles for a Specific Test Case 
With Wall Expansion Angle a ” 90* and Swirl Vane Angle 4 » 45* 
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Figure kl , Probe Interference Effect on Measured 
Velocity Profiles With Wall Expan- 
sion Angle a ■ 90* and Swirl Vane 
Angle ♦ » A5* at Axial Station x/D 
■ 1.0 for (a) Axial Velocity and 
(b) Swirl Velocity 
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9 ) U/Uo 

Figure ^9. Effect of Inlet Reynolds Number on 
Measured Velocity Profiles With 
Wall Expansion Angle a ■ 90** and 
Swirl Vane Angle t ■ kS'‘ t;t Axial 
Station x/0 ■ 2.5 for (a) Axial 
Velocity and (b) Swirl Velocity 
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Figure 50. Measured Velocity Profiles From an Entire 
Diameter Traverse With Wall Expansion 
Angle a ■ 90* and Swirl Vane Angle $ ■ 

A5“ at Axial Station x/D ■ 1.0 for (a) 

Axial Velocity and (b) Swirl Velocity 
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P^ioure 51. Flow F.icillty for Tountcrswl rlf ng and 
Coswirling Fiow for Coannuiar Pipe 
Flow Prediction Test Caoeo (kO) 
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0 0-5 |.0 1-5 2.0 

(a) u/u„ 

Figure f>2. Inlet Velocity Profiles at x/0 ■ 
0.01375 (x - 0.002 m) for the 
Counterswirling Coannular Pipe 
Flow Test Case (a) Axial Velo- 
city and (b) Swirl Velocity 
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Figure 52. (Continued) 
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Figure 5^. Velocity Profiles at x/D ■ 0.27^ (x • A cm) 
for the Counterswirling Coannular Pipe 
Flow Test Case (a) Axial Velocity and (b) 
Swirl Velocity 
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Figure 56. Inlet v'clocity Profiles at x/D ■ 
0.01375 (x - 0.2 cm) for the 
Coswirling Coannul ar Pipe Flow 
Test Case (a) Axial Velocity 
and j) Swirl Velocity 
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Figure 57- 
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Figure 57. (Continued) 
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Figure 58. Velocity Profiles at x/D ■ 0.206 (x »' 3 cm) 
for the Coswirling Coannular Pipe Flow 
Test Case (a) Axial Velocity and (b) Swirl 
Veloci ty 
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C OKCAHCNA STATE UNIVERSITY 

C STILLWATERt OK 74078 

C 

c 

CHAPTER OOOOOOOO PRELIHINARIES OOOOOOOO 
C 

DINENSIIIN HEDM(9) tHEDUNNIf > tHEDNhSI*) •HEDCNU<«> t 
•HFOU< »« .HE0V<9> >HEDU<9> tHE0VT<9) »HCD«JSr(9> > 

«MFPVSr(9> .HEDW8T<9) tMEDVTS«9) .MEDDCC(9)tMEDDET<9> , 

IHEDNNr (9t .HEdHIV<9 ) ,mEDHIP( 9> tHEDnPP(9) tHE0AH(9> 

C 

cohmon 

l/CAL I D/CPI TCH ( 26 > > CDCL r A < 2A ) • CVCLCF ( 36 > 

t/Nf.A8UR/R8CTA<4e>34> •RPNNP8<48>34> tRPCnPU(48t24>»NDATA<48)>NAXJPTf 


» RDNPR0<48) 

*/0E0H/X(4R> >R(34> tXND(48> >RNP<34) •DYP8(24) •0VNP<34> • 
t ‘}NS«34/ tNCTATNtXlNCHS(48> f RINCHS<34) 

t/CALC/ VT0TAL(40t34> >U(4R.34) tVf 48r34)fW(48f34>« 
t VT0TAR<48t34) iU6TAR(48i34t tV8TAR(48>34> •HSTAR<48>34>f 

* PICMrF(4Bf 34>.VELCF«48f34)f DEl.TA(48>34>.RErA(4B>34> t 

»ANONOrl< 48) » IJNEAN(48> tf1ARS(48> tHASFL0(48> tUIN(48> 

REAL HASSfNA'iFLO 
LOGICAL IWRITEtDIAGNS 

C SET lUPITE-.TRUE. FOR WKITINO SOLN. ON DISK STORAGE 

C bet DIAONS-. true. TO ACTIVATE DIAGNOSTIC VRITE 3TATEHENT8 

C-- N3TATN IS NO. OF AXIAL STATIONS FOR WHICH DATA IS 

C SUPPLIED 

C HAXJPT IS NAX. NO. OF RADIAL POSITIONS FOR ALL AXIAL 

C STATIONS 

C CPITCH IS CALIBRATION PITCH COEFF . 

C CDELTA tS CALIBRATION PITCH ANOLFIDCO.) 

C CVCLCF IS CALIBRATION VELOCITY COEFF. 

C NCAL IS NO. OF CALIBRATION POINTS 

C XINCHB IS AXIAL POSITION! IN. > OF EACH PARTICULAR 

C RADIAL traverse 

C .,—NDATA IS NO. OF RADIAL LOCATIONS FOR WHICH DATA IS 

C SUPPLIED FOR EACH PARTICULAR RADIAL TRAVERSE 

C RDNPRS IS MEASURED DYNAMIC PRESSURE ENTERING THE 

C 3WIRLER 

C RINCH3 IS RADIAL POSIT ION! IN. » OF EACH HEASUKEMENT 

C uOCATTON IN EACH PARTICULAR RADIAL TRAVERSE 

C RBETA IS MEASURED FLOW ANnL£<DEG.) WHERE YAW 

C ANGLE BETA - 160.0 - RBCTA 

C RIiNfRS IS MEASURtD VOLTS FOR PNORTH - PSOUTH 

: OIFF. PRESSURE 

C RPCMPW IS MEASURED VOLTS FOR PCENTER - PUEST 

C DIFF. PRESSUKE 

C PICHCF IS REDUCED VALUE FOR PITCH COEFF. 

C delta is reduced value KOR PITCH ANOLE<DEG.) 

C VELCF IS REDUCED VALUE FOR VELOCITY COEFF. 

C BETA IS REDUCED VALUE FOR PROBE YAW ANOLE<DEO.) 

C VTOTAL IS TOTAL VECTOR VELOCITY MAONI TUDE<M/8EC> 

C— — U IS AXIAL VELOCITY 'M/SEO 

C W IS RADIAL VELOCITY IM/SEC) 

C U IS SWIRL VELOCITY (M/SEC) 

C VTSTAR IS DIMENSIONLESS TOTAL VELOCITY 

C USTAR 13 DIMENSIONLESS AXIAL VELOCITY 

; VSTAR IS DIMENSIONLESS RADIAL VELOCITY 

C WSTAR IS DIMENSIONLESS SWIRL VELOCITY 

C ALL PRIMARY USER INPUTS ARE LOCATED HERE 

IWRITE-.TRUE. 

DIAGNS'-.TRUE. 

ALPHA-90. 

PH I -45. ORIGINAL PAGE IS 

OF. POOR QUALITY 


i 


I 


i 


i 

M 

4 


I 


ViriCOM-t .8C-S 
nSTATN-* 
«AXJPT»20 
PArn«'»i.B 


ORIGINAL PA tL 
OF POOR QUALlTy 


tT«48 

JT-J4 

AROC* 1 1. 79 / 0.0*39.37 > 

MGnAt.L«RLAR0C/2.0 

RCAO<?f 20S> HCDM.HCDUnN.HCOUfHCDV.HCnU. 

I HCtiVT.HEtHJST.HCDUSrtHESyST.HE&VTGtHCDDCL.HCOkCT. 

IMP [INM 8 . ME DCHU . HE DHMF • HE fiM t <7 1 ME PH I P > ’ (C DHPP . Hf 0AM 
209 FORMAT <7A4> 

C otNlTIALIZC VARIABLES TO ZERO 

CALL INTT 

C read riVC-MOLE PITOT CALIBRATION DATA 

NCAl *29 

on 10 t>ttNLAL 

Rt A0*9.210) CPITCHl I > .rOELTAC I ) .CVELCFU ) 

10 rriNTINUE 
2tO r0RMAT(3F10.S* 

(F'DIAGNS) URtTE««.400> (CPI TCH< I > . I «1 • 29 > 

IF(OIAONS) aRire<At400> <C0ELTA(I).I«1.29> 

(F<0IAGNS) URITE(4.400> ( COELCF ( I > > I • 1 «29> 

400 FORMATC///. tX. 13(F8. 4. 1X> >//.9X. 12(F8.4) > 

C ....reaO raw measured data to be REOI'CCD 

DO 30 I-l.NSTATN 

RE AD (9. 230) XINCH8( I ) .NDATA( I > tRDNPRS< I > 

JPT5.nPATA« I ) 
no 30 J-l.JPTS 

READ (9. 230) RINCHS( J> iRBETAi I . J> . RPNMPS( I f J> .RPCMPUl I > J) 
20 CONTINUE 
30 CONTINUE 

CONVERT X AND R FROM INCHES TO METERS 

on 39 I-tlfNCTArN 
X I >*XINCHS( I ) *0.0294 
JPT3«NDATA( I > 

DO !2 J-1«JPTS 
Ri )>»RfMCMS( J'*0.0294 

32 C JNTInuE 
19 CONTINUE 
220 F )RMAT(4F10.S) 

230 FORMATdFlO.Sf lIlOtlFlO.9) 

tF<DIAUN6> WRirC(4.470> <NOATA( I > • I-l . NSTATN) 

IF(DIAGNS) URITE(At490> ( X ( I > . I > 1 .NSTATN > 

IF(DIACNS) URITE(4.S00> ( R( J> . J- I . JPTS ) 

DO 37 (-1. NSTATN 

IF(DIAONS) NRITE<A.900> <RBETA( I « J> > J>1 . JPTS) 

IC(DIAONS) URITE<6.900) <RPNMPS< I . J > . J«1 . JPTS ) 

IF(DIAONS) WRITE(A.900> <RPCMPU( I . J ) . J- 1 • JPTS ) 

37 CONTINUE 

)90 F0RMAT<///.40X.l<F8.4.tX)> 

900 F0RMAT(///.20X.10(F8.4) > 

C 

CHAPTER 11111 DATA REDUCTION 111111 
C 

470 FCRMAT(///.40X.1(I9.1X)' 


C CALC PICHCF AND INTERPOLATE FOR DELTA FROM 

C PITOT CALIBRATION CURVE 


IDID'O 

DO 90 I-l. NSTATN 
JPTS-NDATA<I> 

DO 40 J-I.JPTS 

IFiipPCMPU(I.J) .eO. 0.0) .AND. (RPNMPS(I.J) .EQ. 0.0)> 00 TO 38 
F ICHCF ( I . J ' -RPNMPS < I . J ) / < RPCMPU < I . J ) f 1 . E -« ) 

ir<<PICUCF<I.J).0T.3.J99) .OR. < PICHCF < I . J) . L T . -3 . 799 ) > 00 TO 38 
IFUDID .EQ. 0) r'ELTA(I,J)-SPLlNE(CPITCH. 

* CDFLTA.NCAL.PtCHCFd.J)) 

IFdDID .OT. 0) DELTAd, J)-SP(CPITCH.CDELTA. 

» NCAL. PICHCF! I. J) ) 

IDID-1 
GO TO -*0 
38 CONTINUE 
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^hiUlNAL PAGE ta 
OP POOR QUALiTt 


(tCLTAt I • J»«0.0 
WN|TC<*»BSO> I>J 

Ko»H4r<ro<f 'PiCHCF is out of ramoc of cacii»f<»tion at i« 

AND J•'•I3> 

«0 CONTINUL 
^0 'TONTINliC 

C — — INTCRFOLArC FOM UfLCF FROM PITOT CAtlSPATIUM DATA 

IMD -0 

M) SO t«l>MyTArN 
jpr''>>NDATA( I ) 

00 70 J«JtjPT8 

IF(<RPCNPU(IiJ> .EQ. 0.0> .ANO. (RPNNPSdtJ) .CO. 0.0)) 00 TO AS 
IFMASSfOCLTACl.O))) .OT. SO.O) 00 TO JS 
IFdOtO .CO. 9) VCLCFd.J>«3PLIOC(CI)CLTA. 

ICUCLCFtNCALtOCLTAd* J> > 

IFdOIO .OT. 0> '«LCF< l.J)-OP<CtlELTA.C'.»eLCF. 

INCALtDCLTAd. J>) 

I0ID«1 
GO TO 70 
AS rnwTiNuc 

OCLCFd . J)-0.0 
wR1TC<A.890) [>J 

090 FORMA r<20A. 'DELTA IS OUT OF RANGE OF CALlDRATtON OA^A 
tAT I> .I3> AND J«'tl3) 

^0 CONTINUE 
80 CONTINUE 

DO OS l«lfNSTArN 

IFIOIACNS) URITCIA.SOO) (PICHCFd . J) > J-1 . JPT8 ) 

!F<DIACNS) WRITE(AfSOO) < DEL TAd . J) . j-l . JPTS) 

IFIDIAGNS) WRITCIAtSOO' < OELCPd > J) • > JPTS > 

OS CONTINUE 


^ CALC magnitude of total mean velocity vector and 

c U> Vf t W r.OMPONFNTO 


RHO«PATM«( 13 A. 0 / 0 . 102 ) /< 287 . 0 »« TATMf 273 . 0 ) ) 

PI-J. 1 A 1 S 9 

DO too I- 1 .N 8 TATN 

JPTS-NDATAd) 

• 00 90 J-I.JPT 8 
8 ETAd>J)« JAO. -RDETAd.J) 

ir((KFCMPUd>J>.C 0 . 9 . 0 ) .AND. (RFNMPOd > J) .CO. 0 . 0 ) ) RETAd.J)- 0.0 
VTOTALd • J>« 6 QRT(ADS( 2 . 9 /RH 0 »VELCFd. J)«RPi:mPU< !• J)« 133 . 9 ) > 

U< I r J> •VTOTALd > J> • C 03 <DELTAd> J)«PI/ 1 U 0 . 0 ) « 

• C 08 (DETAd. J)«PI/ 180 . 0 ) 

Vdf Jl^VTOTALd* J) • SIN<DELTAd. J)*PI/ 180 . 0 ) 

•J. I. J)*VTOTALd. j) 9 i:OS<DELTA<I. J)*r|/ 180 . 0 ) * 

• 81 N(DC 1 Adt J)«PI/ 180 . 0 ) 

90 i.'ONTINUC 

ICO CONTINUE 

IF<DIAONS) URITC(A>SOOMVTOTALdf J) t J- 1 . JPTS) 

IFIDIAGNS) URITC(A>SOO)(Ud. J) .J>lf JPTS) 

IF<DIA 0 N 8 ) ulRirECA.SOOXVd. J) . J*l. JPT 8 ) 

IF(DIAGNS) WRITC(AtSOO)<Ud.J).J-l>JPTS) 

CHAPTER 222222 AUXILIARY CALCULATIONS 22222 
C 

DP 130 I-l.NSTATN 

C CALC GEOMETRIC OUANTITIES 

JPTS-NDATAI I) 

•JPTSMl-.IPTS-l 
DYPSd **0.0 

DYNP( JPTS>- 2 . 0 «(RLAR 0 E-R( JPT 3 ) ) 

DO no J-l.JPTOMl 
DTNP( J)>R< JAl )-R( J) 

DYP 3 ( JTI >-DYNP< J) 

110 CONTINUE 

DO 115 J-l.JPTS 

SM 8 < J>«O.S«<DYNP( J)+DYP 3 < J) > 

US CONTINUE 

IF<DIAGNS> WRlTE<AtSOO> ( DYNP( J ) . J» 1 . JPT 3 ) 

IF.DIAGNG) URITC(A.SOO) < SNSU ) • J- 1 1 JPTS ) 

FLOU- 0.0 
UHOM- 0.0 
DO 120 J-l.JPTS 
ARDEN-r.HO«R ( J ) tSNS ( J > 
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JHIGINAL P‘ ‘ 
OF POOR QUmLIi V 


WMor<*anon»AKDCN«u( l r ji«b( I • j)«x< jt 
n.OW«l*’LOW»A«cOCN«U( I • J) 

ir<DIAnN«> UKITr(A>900) J*n.OU>ARDCN»UMOniRHO 
VOO njl»«AT<///.I3»4»ri2,4t2)0) 

120 CONTI NIC 

ANGNOm 1 >>wNon 
NA'i3< P»2,O»PI*rU0U 

JHCANI I laAAtSm/IRilOtPIaKLAAOCOJ) 

130 CON ri NIC 

iriDlAONS) WNtrC(4>4S0> (UHeAN(I><{at>NSTATN> 
iriDIAONS) URtT4(4>4S0) (HASSI I ) • !•! •N8TATN) 

C NUNDIrtCNtflONALlZE ^CLOCITICS 

00 130 I-ltNlTATN 
INOI I>«»(I)/<3.0«RLAROC> 

JPTO-NPATAI I » 

U1N( I >«(SnRr<2.0/RH0«R0NPRS( I )«34O,O0) 312/3.44 >«*2 
NASPLOi r )-2.''aPt»r<M0*niN( i >«R3 hallm2/2.o 
00 140 J-1..JPT3 

VrSfARI 1» J>«OTOTAL( I. J)/U1N( t > 

•JETARI 1 f J)«U( I. JI/UINI I > 

ORTARI 1 . J)«VI( 1 . J)/UIN< I > 

WSrARI 1 > J> •U( 1 f J»/U1N< I) 

140 CONTINUli 
ISO CONTINUE 

00 140 J-l.NAXJPT 

RND< J)>R( J)/(2.0«RLAR0E> 

140 CONTINUE 
C 

CHAPTER 3 3 3 3 1 OUTPUT 3 I 1 3 1 J 3 

C 

IM..<IJT. lURITF.l 00 TO 145 
WRirCill) < 

URITEdl) R 
URITEilt) U 
WRITE! ID U 
WRITE! ID W 
143 CONTINUr 

URITE(4.31D 
WRITE! 4 .323) ALPHA 
WMTE!4.330) phi 
WRITE 14.333) R8NALL 
WRI rC!4.340) RLAROE 
WHITe(4.33S> WISC03 
WRirE<4t340) RHO 
C 

CALL WRITE! 1 . 1 .N3TATN. 1 . I T . JT . < .R .AASFLO.HEl'NMF » 

CALL WRITE! I . 1 .NGTATN. 1 . IT , JT . X , R.MASS .HEPN > 

CALL WRITE! 1 . 1 .NSTATN. 1 , IT . JT . X .R. UIN. hEDNIW > 

CALL WRITE! 1 .1 .NSTATN. 1 . 1 T . JT * X .R. UNEAN t HCDUNN ) 

CALL WRITE!! .1. NSTATN. 1 . 1 T . JT > X .R . ANOhON t HEI'AN ) 

CALL PRINT! 1 .1. NSTATN. NAXJPT. IT, JT.X.R.U.HEPU) 

CALL PRINT! 1 .l.NSTATN.NAXJPT.IT. JT.X.R.y.MfDW) 

CALL PRINT! 1. 1 . NSTATN,. NAX. IPT , IT , JT , X ,R .UtHEPU ) 

CAUL PRINT! 1 , I .NSTATN, MAX JPT. IT, JT.X.R, DELTA, HEPPED 
CAl L PRINT! I . 1 .NSTATN. MAX JPT, t T , JT . X ,R. PETA.HEPPET ) 

CALL PRINTM . 1 .NSTATN, MAX JPT . IT . JT . X .R.UTOTAL .HEDUT > 

CALL PRINT! 1 . 1 .NSTATN.MAX JPT . IT . JT . IND.RNP.USTAR . HEPUST ) 

CALL PRINT! 1. 1, NSTATN, MAX JPT. IT, JT.XND.RND.VSTAR.HEPWST) 

CALL PR r NT ! 1 , 1 . NST A TN . MAX JPT . I T , JT . XNP . RND . UST AR . HEDU8T ) 

CALL PR INT! 1 , 1 . NSTATN.MAX JPT • IT . JT . XND .RNP.L'TSTAR .HEPUTS ) 

CALL PRINT (D I, NSTATN, MAX JPT, 1T.J1 . XINCHS.RINCHS .RPNMPS .HEDNMS > 
CALL PRINT! 1 . 1 .NSTATN.MAX JPT. IT. JT.XTNCHS.RINCHS.RPCMPW.HEPCMW) 
call WRITE! 1 . 1 .NSTATN. 1 . IT . jT . XINCHS.RINCHS .RPNPRS .HEDMIP) 

CALL PRINT! 1. 1. NSTATN. MAX JPT. IT. JT, XINCHS.RINCHS, PICHCF.MEDMPP) 
STOP 

FORMAT STArPMENTS 

III rCPMAT! IHI , rST, 'AXI3YMMCTRIC, ISOTHERMAL, OT COMPUSTOR FLOWFIELD 
fMFASI/REMENTS T»3, IJSINO A FIVE -HOLE PITOT PROPF. ' ) 

32C form.'.:!//// , T40, ‘EXPANSION ANGLCCDEC.) • ' , T77, 1PE13 . 3) 

330 FDRMAT!//, T.»0, SWIRL 'JANE ANGLE'PEO.) ^ , T7T , 1PE13 . 3 ) 

373 FORMAT (./, T40 «' INLET RAHIUS(M) « ' , T?’, IPE13 . 3 > 

340 FCf MATl, /,T40« COMPUSTOR RAMUStM) • ' , T77 . I PEI 3 . 3 ) 

=’0RMAr!//,T»0, 'LAMINAR UISCOS ITV ! XG/M/SEC > • ' . T77 , IPE13 . 3 ) 


160 



nnoiinonnofionof- 


» ,0 •//. t40. ' WLNiITY<hO/CU. «> •' .T77. jrcu.i./ V/) 

tlMP 

TUHHOUTImT tNlT 

C«««t«*««a««t««««<v««**«*«««»*««****«*«*««*«*««««***«****««**«*** 

C 

cohhon 

» '*lCA«liK 'A»r r A< *« 1 34 ) t A^NHAS ( 40 1 34 > . AACNAW (40.34). NDAT A ( 40 > . AAX JAT . 
I AriNAA0<40) 

• /OCCA '4 40 ) . A ( 34 > . XND < 40 ) . ANP ( 34 ) . DVA0 « 34 ) . D VNA ( 34 > . 

I b .'U(34).NflTATN.XINCH0<40> .AINCH0(34) 

»/CALC/VTOTAC<40.34>.l)(4a.34>.V<4e.34).W<40.34). 

• VT0TAA(40.34).U0TAA(40.34>.V0TAA(40.34>.USTAA(40.34>. 

I AtCICA(4B.34).V0rLCA(40.34) .D€LTA(40.34).0CTA(40.34>. 

• ANOMUn ( 40 ) . unc AN ( 40 ) . HA08 ( 40 ) . NASA I O(40).IJ|N(4A) 

AC At. 4A‘?0.MA8FIO 

C 

no 30 (ol.NSTATN 
MASFLOt I )«0.0 
MA’7S(t)«0.0 
ArtnNON ( I ) «0 . 0 
iJNCANt T)-0.0 
UtN< I >«0.0 
00 10 J-l.HAXJAT 
yT0TAL(I.J)»0.0 
U< I . Jt«0.0 
V( t. J)-0.0 
U( I . J>«0.0 
VT0TAR(t.J)«O.O 
IJ0TAA< I.J)-0.0 
•^0TAA(I . J)«0.0 
U0TAA( I.J)>0.0 
AACTAd. J)>0.0 
AFTA(I.J)-0.0 
PFNNAAf I . J) *0.0 
AACNAUt I . J)«0.0 
AICHCr(I.J)»0.0 
WCLCAd. J)»0.0 
PCLTAd.J)-0.0 
10 CONTINUC 
30 CONTINUC 
AFTUAN 
FNP 
C 

FUNCTION 8Al.tNC<X. FX. N. 41) 


CUBIC 0AI.INE CUAWC FITTING IN 3 DIMENSIONAL DATA PLANE 
INA'JT UAL DCS I 

X. FX DATA ARRAYS. ONE DIMENSIONAL. X IN INCREASING OAPER 

N NUMBER OF data AOINTS IN 4. MAX 34 

XI POINT OF INTEREST. UHCRE F(X1> IS TO BF FOUND 

RETURN VALUE t 
SALINE OR 3A • F(Xl) 

THIS ROUTINE ACTIVATES ROUTINE ABUILP. M. AND GAUSS. 

FOR INTERPOLATION OF A LARGE NUMBER OF DATA POINTS. FUNCTION 
SPLINE MAY BE CALLED ONLY ONCE . AND .iUBSEOUENT CALLS MAY USE 
ENTRY POINT SP. 


DIMENSION Xd>. FX(1>. A(34.37> 

C CONSTRUCT SPLINE MATRIX 

N1 aNTl 

DO 10 I«l. N 
DO 10 J-1. N1 

>0 Ad.J)*0. 

Ml-N-1 

DO 20 I«2. Ml 

30 CALL ABUlLOtX. FX. A. N. I) 
A(1.1)*H(X.3) 
A(1.3>»-M<X.1)-H<X.2) 
A(1.J>-H(X.1> 

M3-N-3 

A<N.M2)*H(X.M1 > 
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ORIGINAL PACiE ^ 
OF POOR QUALfiY 


.•N.nt >«>H<x>n2> -H(X>ni > 

^•NfNV «M( XtH?) 

C FTNO ‘ICCOND 

iV.i.L liAIJ3'I«A> Hf tm 
•■rITRf 'ir<X. KXt N> <l> 

f iNti ri <n 


1.0 40 t *1 • Al 
n -tfi 

trixt .CQ. ««i>) lio ro so 

tr«Xl .CT. X< 1 > .AND. XI .OT. x<ll)> 00 TO 41 
tr<Xl .OT. X(|> .AND. XI .LT. X(in ) 00 TO 4 t 

40 rONTTNUF 

trfXl .CO. XtN>> 00 TO AO 
aRITE'At 43 > <t 

4 ? FORHAT^' xt»'. 014 . 7 . ' OUT OF INTCRF 01 .ATI 0 N RANGE. RETURNFP VALUE 

»•«)'' 

'TPI.IME-). 

sror 

41 i;ONTTNIIE 

n«ui 

Nt*H<X. t> 

HX»*« tl >-<t 
HX ’•XI -X( I > 
f*< l/Mf -NI«HX 

r <l*F <t «A( I .Nt ) 

■irO.HX - MI»HX? 

' <i •' Fvt .srotAi n ,Ni > >/A. 

.PL INE •' FX< t **HXfr << 1 1 ) «HX?> /Ht4FX t 
IP^OPI.INE 
.•'ETUPN 
r 

SC C ON T tout 

npi.iNt«rx' t ) 

"P'SPt INF 
HFTIIRN 
r 

CONTINUE 

• 3 PI.INC»FX<N> 

IP .SPLINE 

RETURN 

CND 


FL'NCTION H< X. I ) 


0 calculate nr. T. x which is usijailt called .^s h. 

riMENSION <<l) 

n-i»i 

H«X< I I ) -X« I > 

return 

END 

jUBBOUTINF ABUILD(X. f. a. n. i» 


C CONSTRUCT SPLINE MATRIX FOR FINDING CND OERIVATIUES. 
ril**»'***t**** 1 l*M*f**»t *»$****«*»*$«*** 

DIMENSION X<t». Fil). A(2A.r7) 

lMI-t-1 

Il-IM 

Nl.N+l 

S’’0»H(X. I ) 

HIM1»H»X.IM1> 

A< I. IMl >«HIMI 
.T<l.H«2.4(HlMt+3T0> 

A(I.m«3T0 

A<I.N1)-< <F<Il)-F(I>>/3Tn - (r(I)-F<IMl))/HIHI )t4. 

return 

END 


jUPROUriNE '5AU33<A. X, M ) 

gausf-joroan elimination 
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n r* 


C*»4» *«•«*» »«*****t*t***«*<*** 

DtfiCNniON A<2i>27' 

DO 1 L«li K1 
Ul-LM 

00 T 1-Ll. K 
C0NST<A( ItL)/A\Lri.) 

DO 3 J«L» M 

3 A* tf J)-C0N0T*A<L* J> 

DO A t«l> M 
n«i^t 

DO A L«ll> HI 
rONST«A<I»L>/A<L.L> 

DO A J«I. H 

A An>J)-*A< If J>>CONST«A<L»J> 

00 10 I-l> K 
A<IiH>>A<I>n)/A(I>I) 

10 A(I>I»«1. 

RETURN 

END 

C 

SUDROUTINE PRINT ( ISTART • J8TART »NI >NJ • I T t JT t X t V t PHI r HEAD > 


DIHENSION THK ITt JT) tX< IT > t T < JT ) * HEAD< 9 ) 1 3T0RE < A8 > 

ISMP-t 

J8KIP-1 
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